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Resumen 



En la presente tesis se reporta un conjunto de invest igaciones teoricas realizadas por el autor 
en el periodo de 1991 a 2005, es decir, con posterioridad a la defensa de su Tesis de Candidato 
a Doctor en Ciencias Fisico - Matematicas en 1990 en el Instituto de Fisica P.N. Lebedev de 
Moscu. Los sistemas fisicos que se estudian abarcan modelos de quarks de bariones, clusters 
atomicos, atomos confinados en trampas, asi como electrones y excitones en puntos cuanticos. 
El comun denominador de todos ellos es que son descritos por la Mecanica Cuantica no rel- 
ativista y que se componen de tres o mas particulas. Incluso para el sistema mas pequeno 
considerado, el de tres particulas, no existen soluciones analiticas exactas que permitan hacer 
comparaciones o den una vision cualitativa de sus propiedades. Por otro lado, en un prob- 
lema donde intervienen decenas o centenares de particulas aparecen fenomenos nuevos como 
excitaciones colectivas ( "plasmones" ) y otros que son determinantes en la explication de sus car- 
acteristicas. Los metodos utilizados en el estudio de estos sistemas van desde la diagonalizacion 
directa acompanada por el algoritmo de Lanczos y algunos metodos analiticos (la expansion en 
potencias de 1/D, la cuantizacion semiclasica en conjuncion con transformaciones canonicas, 
etc) en los sistemas mas pequenos (menos de 6 particulas), hasta metodos aproximados en los 
mas grandes (de 7 a cientos de particulas) como, por ejemplo, los aproximantes dobles de Pade, 
el metodo variacional de Monte Carlo, Hartree-Fock y la aproximacion de fase aleatoria (RPA) 
en su version de Fisica Nuclear para sistemas finitos, la aproximacion BCS y el algoritmo de 
Bethe-Goldstone de la Fisica Nuclear. El rasgo comun de todos estos metodos es su caracter no 
perturbativo y la necesidad de recursos computacionales relativamente importantes. El autor y 
sus colaboradores han utilizado diversas estrategias para optimizar los procedimientos reducien- 
do asi el costo computacional. Un ejemplo ilustrativo es el caso estudiado de 12 electrones en un 
campo magnetico intenso donde la matriz hamiltoniana, de dimension superior a 170 millones, 
es imposible de diagonalizar directamente y, sin embargo, sus autovalores mas bajos fueron 
hallados utilizando una combination del algoritmo de Lanczos con la teoria de perturbaciones 
degenerada de segundo orden. Las propiedades fisicas estudiadas abarcan los espectros de en- 
ergias, las densidades, la textura del espin, el radio cuadratico medio, la absorcion y emision 
de luz en distintas frecuencias, la dispersion inelastica de la luz y la evolution dinamica hasta 
llegar a un regimen estacionario. En algunos de los casos, el autor y sus colaboradores han sido 
los primeros, o incluso los unicos, que han descrito determinadas caracteristicas de estos sis- 
temas tan complejos. Para citar un ejemplo, no existen en la literatura cientifica otros calculos 
de dispersion inelastica (Raman) de luz en puntos cuanticos semiconductores con mas de 40 
electrones, aiin cuando los experimentos reportados estiman que el niimero de electrones puede 
llegar hasta 200. En general, la novedad del trabajo que se presenta consiste en la aplicacion por 
primera vez de metodos, desarrollados con anterioridad o utilizados en otras ramas de la Fisica, 
al estudio de problemas en que intervienen desde tres hasta decenas o cientos de particulas y 
cuya fenomenologia es nueva y esta en plena discusion en la literatura cientifica del momento. 
Utilizando estos metodos se han hecho predicciones importantes sobre propiedades de los sis- 
temas fisicos estudiados, que pueden ser sin dudas verificadas experimentalmente. Entre ellas 
podemos citar, por ejemplo, la existencia de resonancias gigantes en la absorcion optica de 
puntos cuanticos fuertemente excitados o la aparicion de discontinuidades en la densidad de 
estados de puntos cuanticos con pocos electrones bajo campos magneticos intensos. 
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Capitulo 1 
Introduction 



1.1 Unidad de la tesis 

La presente tesis de doctorado en ciencias ha sido disenada a posteriori, es decir despues que 
los trabajos que la constituyen fueron concluidos. A lo largo de estos 15 anos el autor ha 
abordado problemas con fenomenologia diversa que abarcan sistemas de la Fisica subnuclear, los 
conglomerados (clusters) de atomos, el condensado de Bose, los denominados atomos artificiales 
(electrones en un punto cuantico) y sistemas que exhiben una gran analogia con los nucleos 
atomicos (multiexcitones en puntos cuanticos). Fue una sorpresa para el autor descubrir la 
unidad en esta multitud y darse cuenta que la misma se basa en el impulso recibido durante el 
trabajo de candidatura en Moscu en el intervalo de 1986 a 1990. 

Asi, en la candidatura que se culmino en 1990, siguiendo una sugerencia del Prof. D.A. Kirzh- 
nits, se abordaron principalmente problemas de 3 y 4 particulas en la Fisica Nuclear y Atomica. 
El foco de atencion eran los autovalores mas bajos de energia y, muy especialmente, la obten- 
cion de relaciones universales, es decir independientes del potencial de interaction entre las 
particulas. 

En el trabajo posterior, que constituye el grueso de la presente tesis, se amplio el rango de los 
sistemas considerados, llegando a estudiar problemas con 400 particulas, se diversificaron los 
metodos utilizados en este estudio y se amplio sobremanera la fenomenologia tratada, que ahora 
comprendio no solo otras propiedades estaticas de estos sistemas sino, ademas, algunos aspectos 
de la optica lineal y no lineal de los mismos y de la dinamica en condiciones de transiente. El 
espiritu, sin embargo, se puede afirmar que es el mismo: el estudio de sistemas de 3 o mas 
particulas utilizando metodos esencialmente no perturbativos. 

El co nj unto de articulos reportados como soporte de la tesis [I] - [HH| responden a esta unidad. 
Han quedado excluidos solo 3 trabajos, concluidos tambien en este periodo, que abordan, 
respectivamente, la distribution de corrientes en el efecto Hall cuantico entero |40j . la analogia 
entre un problema de Dirichlet no lineal y las ecuaciones de Newton con friction |3T] y un 
comentario sobre las trampas atomicas [121 ■ El autor considera que estos 3 articulos se apartan 
de la linea general. 
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1.2 Complejidad de la tematica abordada 

El problema (clasico o cuantico) de una particula moviendose en un campo exterior tiene 
grandes atractivos. Por un lado, es el sistema conceptualmente mas simple, donde muchas de 
las regularidades de la Fisica se ponen de manifiesto. El analisis cualitativo de muchos pro- 
cesos fisicos se basa en un cuadro de particulas independientes, donde las interacciones entre 
los cuerpos se desprecian en una primera aproximacion. Asi, por ejemplo, muchas propiedades 
electronicas de los semiconductores pueden comprenderse sobre la base de este esquema. Por 
otro lado, se conoce una multitud de soluciones analiticas correspondientes a distintos poten- 
ciales externos con cierta simetria. Estas soluciones analiticas ofrecen un punto de partida para 
analisis perturbativos o sirven de bases para la description numerica de problemas en que no 
existe la solution exacta. El caso de dos cuerpos interactuando a traves de fuerzas centrales se 
reduce al de una particula en un campo exterior. 

En contraposition, los problemas fisicos en que participan mas de dos cuerpos interactuantes 
son reconocidos como extremadamente complejos. Por ejemplo, el problema clasico gravitatorio 
de tres cuerpos: sol - tierra - luna, con una historia cientifica de cientos de anos, permanece 
retando a los cientificos hasta nuestros dias jlSj y, aiin cuando las ecuaciones de Newton para 
tres particulas se pueden integrar numericamente, la comprension cualitativa de su dinamica 
dista mucho de ser completa. 

En el ambito cuantico, a la complejidad natural del problema de tres cuerpos se le suma el 
hecho de que a veces la Mecanica Cuantica contradice la intuition clasica. Un ejemplo muy 
ilustrativo es el denominado efecto Efimov 03]. Para dos particulas con interaccion de corto 
alcance atractiva que escribiremos g V(r), donde se ha separado explicitamente la constante 
adimensional g, existe un valor umbral go en si cual aparece un estado ligado. Es decir cuando 
g < go no hay estado ligado, mientras que cuando g > go si lo hay. Consideremos ahora el 
sistema de tres cuerpos con interaccion por pares cuando g toma el valor umbral. En ese caso 
existe una sucesion infinita de estados ligados cuyas energias tienden a cero. La mayoria de 
esos estados desaparece al aumentar g, quedando solo el base. Esto es algo dificil de asimilar 
desde el punto de vista intuitive 

Grandes dificultades aparecieron en los anos 60 durante el estudio de los estados de tres cuerpos 
con energia mayor que cero (estados dispersivos) vinculadas al comportamiento asintotico de 
la funcion de onda. La motivation en este caso fueron los experimentos de dispersion neutron - 
deuterio en Fisica Nuclear o electron - atomo de Hidrogeno en Fisica Atomica. Estas dificultades 
fueron resueltas con la reformulation de la ecuacion de Schrodinger como un sistema de tres 
ecuaciones acopladas que se conocen como ecuaciones de Fadeev La demostracion de 

que se puede construir una base completa de funciones para el problema de dispersion fue 
posible en este marco. La generalization al problema de N cuerpos interactuantes se hizo casi 
inmediatamente despues 

Con respecto a la existencia de soluciones analiticas exactas en un sistema cuantico de N 
particulas debemos decir que se conocen un conjunto de soluciones en el caso de movimiento en 
una dimension [46^ , las cuales se relacionan con propiedades de simetria de los correspondientes 
hamiltonianos. Para el caso de movimiento en dimension mayor o igual a dos practicamente no 
se conocen soluciones exactas, excepto el caso (trivial) de interaccion con potenciales armonicos, 
el cual con una transformation de coordenadas se reduce al problema de N osciladores inde- 
pendientes. 
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Por otro lado, en la solution numerica del problema de N = 1 particula, el procedimiento 
usual es construir una base de funciones y expresar el hamiltoniano del problema como una 
matriz en esta base. El mimero de funciones, Ni part , coincide con la dimension de la matriz 
hamiltoniana. Usualmente los primeros autovalores y autovectores de esta matriz, digamos los 
20 primeros, son responsables de muchas de las propiedades fisicas. Para hallar con suficiente 
precision estos 20 autovalores uno toma una base con Ni part = 200, por ejemplo. La dificultad 
que surge al utilizar este metodo en el caso de N > 2 particulas es de naturaleza puramente 
combinatoria. Con Ni part estados uniparticulares uno puede formar Ni part \ / ((Ni part — N)\N\) 
estados del sistema de N particulas. En el caso Ni part = 200, por ejemplo, el niimero de 
estados es mayor que 64 millones solo en el caso N = 4. Y la diagonalizacion de una matriz 
de este tamano es practicamente imposible. Uno puede reducir la dimension utilizando leyes 
de conservation, como la del momento angular, que nos permiten trabajar en subespacios mas 
pequenos. Pero, en la practica, los sistemas con 6-8 particulas son casos limites. 

La mayoria de los problemas tratados en la presente tesis clasifican en el grupo donde la 
diagonalizacion exacta no es posible. Por esa razon, los metodos utilizados hacen determinadas 
aproximaciones. Sin embargo, en todos los casos el tratamiento de la interaction entre las 
particulas es exacto en todo orden de aproximacion, es decir son metodos esencialmente no 
perturbativos. 

1.3 Novedad del trabajo de tesis 

En general, la novedad del trabajo que se presenta consiste en la aplicacion por vez primera 
de metodos, desarrollados con anterioridad o utilizados en otras ramas de la Fisica, al estudio 
de problemas en que intervienen varias particulas y cuya fenomenologia es nueva y esta en 
plena discusion en la literatura cientifica del momento. Utilizando estos metodos se han hecho 
predicciones importantes sobre diversas propiedades de los sistemas fisicos estudiados, muchas 
de las cuales pueden ser, sin dudas, verificadas experiment almente. 

Asi por ejemplo, en la Ref. pOj se estudian sistemas de varias particulas en una dimension. El 
metodo aqui es la supersimetria j!7], que en las Teorias de Campos tiene un significado mas 
profundo, pero en la Mecanica Cuantica puede verse simplemente como una metodica para 
construir hamiltonianos casi isoespectrales. El caso de una particula en un potencial externo 
se habia tratado en el trabajo 0S|- En pQ se construyeron potenciales de interaction entre N 
particulas en una dimension para los cuales se conoce exactamente la energia y la funcion de 
onda del estado base. 

La supersimetria en conjuncion con la cuantizacion semiclasica y las transformaciones canonicas 
se utilizo en los articulos 001313 para obtener los estados excitados de moleculas lineales y los 
estados de tres particulas con momento angular total igual a cero. Para una particula en una 
dimension ya se habia reportado que una regla de Bohr-Sommerfeld modificada a sugerencia de 
la supersimetria era exacta en cierta clase de potenciales jlHl EH EIJ • Pero la extension al caso 
de varias particulas de la receta de cuantizacion semiclasica modificada por la supersimetria 
y la utilization de las transformaciones canonicas para reducir el hamiltoniano a una forma 
normal, facilmente cuantizable, fue realizado, por primera vez, por nosotros. 

Los trabajos E], jlj y [3 EH 13 USj son reminiscentes directos del trabajo de candidatu- 
ra. En ellos se empleo el denominado metodo 1/D consistente en utilizar el inverso de la 
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dimension espacial (o del momento angular) como parametro no perturbativo en la ecuacion 
de Schrodinger p2j . En el articulo [2] se relacionaron los tamanos de los hadrones (compuestos 
de tres quarks, entre ellos el proton, por ej.) y de los mesones (compuestos de un quark y un 
antiquark) y se hizo una comparacion cualitativa con resultados experimentales. Las relaciones 
se obtienen en la aproximacion principal del metodo 1/D. En la Ref. [3] el metodo se uti- 
lizo para obtener el espectro de mesones en un modelo relativista con interacciones armonicas, 
propuesto por Moshinsky. En el trabajo [I] la energia de ligadura y el espectro de excitaciones 
de conglomerados de atomos de gases nobles que contienen hast a 13 atomos fueron obtenidas 
y comparadas con simulaciones numericas que utilizaron el mismo potencial de Lennard- Jones. 
En los articulos EH EH E3 el interes se desplazo hacia sistemas del estado solido como 
el exciton en un pozo cuantico P , pocos electrones en un punto cuantico [TUl EH y anyones 
confinados [15 . Se puede decir que el aporte principal en todos estos casos es la utilizacion 
del cuadro cualitativo que brinda el metodo 1/D al analisis de propiedades fisicas tales como 
el espectro energetico, el radio cuadratico medio y otras de sistemas que contienen hasta 13 
particulas. La aplicacion de este metodo a problemas con tres o mas particulas fue iniciada por 
Herschbach en el estudio de iones de tipo He [53] 152] . Mas adelante se extendio a los atomos 
multielectronicos [HI]. Pero la utilizacion en otros sistemas de la Fisica Atomica, subnuclear y 
de la Fisica del Estado Solido es de nuestra autoria. 

El metodo de los aproximantes dobles de Pade fue utilizado por primera vez en los traba- 
jos [TTJ H21 El HEl HZ| para obtener aproximaciones analiticas a la energia del estado base y 
los primeros estados excitados de electrones en puntos cuanticos y atomos con espin entero 
(bosones) confinados en trampas. Los aproximantes de Pade fueron ampliamente empleados en 
otras areas de la Fisica como, por ejemplo la teoria de transiciones de fase. Nosotros lo conoci- 
mos a traves del trabajo de MacDonald y Ritchie [HE] donde se estudian los niveles del exciton. 
La idea es sencilla: si se conocen las series perturbativas de la energia alrededor de los puntos 
Po y Pi, entonces con el aproximante se construye una aproximacion analitica a la energia en 
todo el intervalo (Po,Pi). En el caso de electrones en un punto cuantico, por ejemplo, donde 
P es la fortaleza de la interaction de Coulomb entre electrones, el caso P — > se puede tratar 
por teoria de perturbaciones y el caso P — > oo tambien se puede resolver semi-analiticamente 
y corresponde a una version finita del llamado cristal de Wigner. En el articulo [IT] , sistemas 
de hasta 5 electrones en un punto cuantico fueron tratados por este metodo. Mas tarde el 
numero de electrones se elevo hasta 210 [Hj. Por otro lado, en las Refs. [13 E] se estudiaron 
sistemas de bosones. En el trabajo [TE] se hizo uso de los dos limites analiticos de la ecuacion de 
Gross-Pitaevskii: la teoria de perturbaciones cuando la interaction entre atomos es muy debil 
y el denominado limite de Thomas- Fermi cuando la interaction es muy fuerte [55] . Con ellos 
se construyo el aproximante de Pade. En la Ref. [T7] se estudio la energia del estado base de 
hasta 200 bosones confinados en una trampa y con interaction de Coulomb entre ellos (iones). 
Como comparacion, se hizo un calculo por Monte Carlo variational [HZj de la misma magnitud 
ofreciendo muy buenos resultados. Por ultimo, en el trabajo [12] se hizo uso de dos desarrollos 
en serie para construir una teoria de perturbaciones "renormalizada" que en espiritu es similar 
a los Pade. El sistema estudiado fue la impureza donora en un pozo cuantico de barrera fini- 
ta. En todos estos casos, los aproximantes de Pade ofrecieron aproximaciones analiticas a las 
propiedades exactas de los sistemas fisicos con errores que no rebasaron el 3 %. 

En los trabajos [TH] EH 123 E3 EE] se hace uso extensivo de la denominada funcion BCS 
para describir el apareamiento de fermiones. Esta funcion es la base de la teoria clasica de 
la superconductividad y ha sido tambien muy empleada en la teoria del micleo atomico [58J. 
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Nosotros la utilizamos por primera vez para estudiar sistemas con un numero finito de pares 
electron-hueco confinados en un punto cuantico. Los mismos se logran bombardeando el punto 
cuantico con laseres de gran intensidad. En el articulo ^H] el enfasis se puso en la region 
de alta densidad de pares y la transition a un regimen del tipo BCS cuando la densidad 
disminuye. El sistema mas grande analizado contenia 210 pares, es decir mas de 400 particulas. 
En la Ref. |2Sj estudiamos la textura del espin y la luminiscencia cuando el punto cuantico se 
halla bajo la influencia de un campo magnetico fuerte. Debido al apareamiento y el caracter 
cuasi-bosonico de los pares, la luminiscencia manifiesta un reforzamiento coherente. En el 
trabajo [35] estudiamos los efectos de la interaction de Coulomb sobre la position de las lineas 
de luminiscencia coherente. En el articulo JSHj se estudio la dinamica de un punto cuantico 
acoplado resonantemente a una microcavidad optica. Experimentos en sistemas de este tipo 
se han realizado recientemente buscando obtener fuentes que emitan pulsos con un solo foton 
|HH1 EDI EI] • Por ultimo, en el trabajo [IB] se propuso un metodo estocastico para corregir la 
deficiencia de la funcion BCS en sistemas finitos, es decir el hecho de que no conserva el numero 
de particulas. 

En la funcion BCS el apareamiento es solo entre pares electron-hueco y se pueden describir solo 
sistemas neutros, en que existe la misma cantidad de electrones que de huecos. La complejidad 
de este esquema es analoga a la de las ecuaciones de Hartree-Fock, lo cual recuerda que BCS 
es tambien una teoria de campo medio. Una aproximacion un poco mas consistente, que tiene 
en cuenta todas las correlaciones de pares, se basa en las denominadas ecuaciones de Bethe- 
Goldstone, ampliamente utilizadas en la Teoria de Nucleos [H2]- En el trabajo [20J se hace, por 
primera vez, uso de estas ecuaciones para estudiar sistemas de hasta 12 pares y sistemas no 
neutros (hasta 6 electrones y 2 huecos). La complejidad numerica no permitio abordar sistemas 
mas grandes. 

En los trabajos [2111211123 123 IS] se em pl eo 1 & diagonalizacion exacta del hamiltoniano cuantico 
para obtener diversas propiedades de sistemas relativamente pequenos. En [21] se calculo la 
absorcion y luminiscencia interbandas de puntos que contienen hasta 3 electrones y en presencia 
de campos magneticos muy fuertes. Estos campos permiten obtener convergencia numerica 
con una base relativamente reducida de funciones multielectronicas (de dimension < 3000). 
En el trabajo [21] se estudio la conductancia de un punto cuantico que contiene un maximo 
de 6 electrones. Experimentos de este tipo fueron realizados recientemente jH3]- El calculo 
independiente de la conductancia y de la densidad de niveles del punto permitio concluir que 
midiendo experimentalmente la primera se podria estimar la segunda. El tamano de las matrices 
se extendio hasta alrededor de 40,000 por lo que fue preciso recurrir a algoritmos del tipo 
Lanczos para obtener los autovalores mas bajos. En la Ref. [26J la diagonalizacion exacta 
combinada con Lanczos se utilizo para obtener la absorcion intrabanda del biexciton en un 
punto cuantico. En el articulo [2H] parametrizamos la densidad de niveles excitados en un 
punto con 6 electrones hasta energias de excitation del orden de 1 meV en campos magneticos 
de alrededor de 8 T. Las matrices diagonalizadas alcanzaron dimension de 850,000, las cuales 
constituyen nuestro record personal. En el trabajo [3T] calculamos la absorcion intrabanda de 
excitones multiplemente cargados, hasta el sistema con 5 electrones y un hueco. 

Los trabajos [221123 estan dedicados a la absorcion intrabanda en sistemas de decenas o cientos 
de pares electron-hueco en un punto cuantico, mientras que los articulos [23 EDI E21 EH EH 
EH] abordan la dispersion inelastica de luz (efecto Raman) en puntos cuanticos. El comun 
denominador de todos ellos es que utilizan aproximaciones de campo medio tanto para el 
estado base (Hartree-Fock) como para los estados excitados (la denominada aproximacion de 
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fase aleatoria, RPA). En dichos articulos se utiliza por primera vez en sistemas de excitones la 
version de la RPA para sistemas finitos empleada en micleos |SH1- En las Refs. |2"2"]l2"3~] se muestra 
que los sistemas multiexcitonicos son analogos a los micleos atomicos y exhiben "resonancias 
gigantes", es decir estados que concentran practicamente toda la absorcion. Experimentos en 
esta direction no existen, pero tampoco serian dificiles de realizar. Se requeriria un bombeo 
fuerte para producir gran cantidad de pares en el punto (la energia del laser debe ser mayor que 
la brecha del semiconductor) y monitorear la absorcion intrabanda (es decir en el intervalo de 
energias del orden de los 10 meV). En el articulo [21] tambien se estudian sistemas neutros con 
decenas de pares pero se calcula la dispersion inelastica de luz. Tampoco existen ni experimentos 
ni calculos parecidos a los nuestros. En las Refs. jSOJ [321 ES E3 EH] se calcula el efecto Raman 
en puntos con decenas de electrones. La motivation son varios experimentos de los anos 90 
en los que se mide la dispersion Raman en puntos con hasta 200 electrones para detectar 
excitaciones colectivas de carga y espin en estos sistemas. Resaltemos que en la literatura 
cientifica el sistema mas grande para el que se ha calculado la section de dispersion Raman (con 
determinadas limitaciones en el algoritmo de Hartree-Fock, etc) es el punto con 12 electrones 

ESI- 

Finalmente, en el trabajo se abordan problemas cuya matriz hamiltoniana es, en principio, 
extragrande, de dimension 170 millones. Estas matrices surgen en la descripcion variacional de 
puntos cuanticos con 12 electrones. Se propone un algoritmo basado en la combination de la 
teoria de perturbaciones a segundo orden con el metodo de Lanczos para evitar la diagonal- 
izacion de matrices tan grandes y obtener el espectro de excitaciones a bajas energias. 

1.4 Estructura de la tesis 

La tesis se compone de un capitulo introductorio, 8 capitulos de resultados y un capitulo final 
donde se ofrecen las conclusiones y las perspectivas hacia el futuro. 

Los resultados se han agrupado de acuerdo al metodo utilizado en la solution del problema 
fisico mas que por la fenomenologia abordada. 

En el Capitulo II (Sec. 12.11 del resumen) se reportan resultados de 7 articulos que tienen en 
comun la utilization del metodo 1/D para estimar el espectro de energias y otras propiedades es- 
encialmente no perturbativas de sistemas cuanticos relativamente pequenos (hasta 13 particulas). 
Se aborda el estudio del radio de los hadrones, la energia de ligadura y el espectro de excitaciones 
de conglomerados (clusters) atomicos y de puntos cuanticos con varios electrones. 

En el Capitulo III (Sec. 12.21 del resumen) se incluyen resultados que tienen que ver con la 
supersimetria en la Mecanica Cuantica y la receta de cuantizacion semiclasica mejorada por la 
supersimetria. Son resenados 5 articulos que tienen que ver con esta tematica, en particular con 
la obtencion del espectro de excitaciones de moleculas triatomicas en una y tres dimensiones. 

En el Capitulo IV (Sec. 12.31 del resumen) se reportan resultados de 4 articulos que utilizan los 
aproximantes dobles de Pade para estimar la energia de puntos grandes (cientos de electrones) 
o de trampas atomicas, donde son confmados desde decenas de miles hasta millones de atomos. 

En el Capitulo V (Sec. 12.41 del resumen) se agrupan resultados que tienen que ver con el uso 
de la funcion BCS y la ecuacion de Bethe-Goldstone. Ambos metodos describen de forma 
aproximada el apareamiento de fermiones. Son utilizados en la descripcion de las propiedades 
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opticas y electronicas de sistemas de excitones en puntos cuanticos semiconductores. Al capitulo 
tributan 5 articulos. 

El Capitulo VI (Sec. 12.51 del resumen) trata sobre la utilization del metodo de Monte Carlo, 
del cual son ejemplos los problemas abordados en los dos articulos que se resenan, es decir, la 
obtencion de la energia de un sistema con cientos de particulas y la estimation de la energia de 
apareamiento en nucleos.. 

El Capitulo VII (Sec. 12.61 del resumen) reune 5 articulos que aplican la diagonalizacion exacta 
y el algoritmo de Lanczos. La fenomenologia aqui es amplia y tiene que ver con la optica 
lineal de puntos cuanticos, los procesos de conduction (a traves de tunelamiento resonante) y 
la pametrizacion de la densidad de estados cuanticos en estos sistemas. 

El Capitulo VIII (Sec. 12.71 del resumen) reune 8 articulos que utilizan metodos del tipo Aprox- 
imacion de fase aleatoria (RPA) para obtener los estados excitados de sistemas relativamente 
grandes. La fenomenologia en este caso esta relacionada con los estados colectivos que se 
distinguen en los procesos de absorcion optica y con la dispersion inelastica de luz. 

Por ultimo, en el Capitulo IX (Sec. 12.81 del resumen) se reseria un articulo que hace uso de la 
combination de la teoria de perturbaciones degenerada con el metodo de Lanczos para abordar 
problemas donde la matriz hamiltoniana es de dimension extragrande. 

En cada capitulo se anexan uno o varios de los trabajos mas representatives publicados por el 
autor con el objetivo de ampliar y detallar el resumen contenido en el mismo. En total se han 
anexado 23 articulos. 

En las Conclusiones se listan los que, a juicio del autor, son los principales resultados de la tesis 
y se senalan las perspectivas inmediatas de continuar el trabajo. 

La bibliografia contiene los 39 articulos del autor que constituyen el cuerpo de la tesis. Ademas, 
se citan 7 de los restantes trabajos del autor y 74 referencias basicas sobre los metodos princi- 
pales utilizados en la tesis y algunos articulos experimentales de gran relevancia. Si incluimos 
las citas bibliograficas de los anexos, el numero de citas supera las 300. 

1.5 A vales de la tesis 

El principal aval de la tesis es la publication en el periodo de 1991 a 2005 de 37 articulos en 
revistas de impacto y amplia circulation internacional [I] [211, [SO] E2|, [Mj GH] y 2 articulos 
en la Revista Cubana de Fisica [23 EH] • Ademas, se han dirigido 3 Tesis de Maestria y 2 de 
Doctorado en Fisica. Una tercera Tesis de Doctorado ya concluida se defendera proximamente 
y una cuarta esta en ejecucion. Los resultados se han presentado en conferencias cientificas 
nacionales e internacionales, han constituido el cuerpo de proyectos de investigation dirigidos 
en Cuba o en el extranjero y fueron objeto de reconocimientos como, por ejemplo, premios 
anuales de la Agenda de Energia Nuclear en 1999 y 2003 y premios nacionales de la Academia 
de Ciencias de Cuba en 1993, 2001 y 2002. 

1.5.1 Participacion en Conferencias 

• Simposio de la Sociedad Cubana de Fisica, La Habana, 1991 
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• Workshop on Condensed Matter, Atomic and Molecular Physics, Trieste, 1991 

• Workshop on Condensed Matter, Atomic and Molecular Physics, Trieste, 1992 

• Oaxtepec Conference on Nuclear Physics, Oaxtepec, Mexico, 1993 

• Course on Geometric Phases, ICTP, Trieste, 1993. 

• Encuentro de Geometria Diferencial en Fisica, Universidad de Los Andes, Bogota, 1994 

• Escuela de Fisica Teorica, Universidad de Antioquia, Medellin, 1994. 

• Adriatico Research Conference on Chaos in Atoms and Molecules, ICTP, Trieste, 1995. 

• Escuela de Fisica Teorica, Universidad Pedagogica, Bogota, 1995. 

• Escuela de Fisica de la Materia Condensada, Univ. de los Andes, Bogota, 1996. 

• Conference on few-body systems near criticality, European Center for Theory (ECT), 
Trento, 1997. 

• Conferencia Internacional "Cibernetica, Matematica y Fisica 97" . La Habana, 1997. 

• Conference on Perspectives in Hadronic Physics, ICTP, Trieste, 1997. 

• Conferencia Internacional "Cibernetica, Matematica y Fisica 99". La Habana, 1999. 

• Workshop on Condensed Matter, Atomic and Molecular Physics, ICTP, Trieste, 1999. 

• Escuela de Fisica de la Materia Condensada, Bucaramanga, 2000. 

• Tercer Taller Caribeho de Mecanica Cuantica, Particulas y Campos, La Habana, 2000. 

• Simposio de la Sociedad Cubana de Fisica, La Habana, 2002. 

• Workshop on mesoscopic systems and Coulomb interactions, ICTP, Trieste, 2002. 

• Third Stig Lundqvist Conference on Advancing Frontiers of Condensed Matter Physics, 
ICTP, Trieste, 2003. 

• 6ta Escuela de Fisica de la Materia Condensada, Medellin, 2004. 

• 17th Simposio Latinoamericano de Fisica del Estado Solido, La Habana, 2004. 

• X Simposio de la Sociedad Cubana de Fisica, La Habana, 2005. 
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1.5.2 Trabajos de Diploma, M.C. y Dr. en Fisica relacionados con 
la tesis 

• Rene Martinez, "Radio de los hadrones debidos a la interaction fuerte a partir del metodo 
1/d" (Tesis de Diploma, UH, 1993). 

• David Leal, "Energias de enlace y excitation en clusters de Lennard- Jones a partir del 
metodo 1/d" (Tesis de Diploma, UH, 1993). 

• Ricardo Perez, "Dos anyones en un potencial de Coulomb" (Tesis de Diploma, ISCTN, 
1994). 

• Boris Rodriguez, "Sistemas pianos de tres cuerpos con interacciones de Calogero y un 
campo magnetico" (Tesis de M. C, UNALMED, 1997). 

• Ricardo Perez, "Sistemas de pocos anyones en un punto cuantico" (Tesis de M. C, ISCTN, 
1998). 

• Alain Delgado, "Resonancias dipolares gigantes en puntos cuanticos" (Tesis de M. C, 
ISCTN, 2000). 

• Boris Rodriguez, "Nucleos artificiales: un estudio teorico sobre las correlaciones elec- 
tron-hueco para el estado basico de un sistema finito confinado en un punto parabolico 
bidimensional" (Tesis de Dr. en Fisica, UdeA, 2002). 

• Ricardo Perez, "Espectros de energias y transiciones internas en sistemas multiexcitonicos 
confinados" (Tesis de Dr. en Fisica, defendida en la UdeA, 2003). 

• Alain Delgado, "Dispersion inelastica de luz por excitaciones electronicas en atomos ar- 
tificiales de GaAs" (Tesis de Dr. en Fisica, proxima a defenderse, ISCTN). 

• Herbert Vinck-Posada, "Optica no lineal de un punto cuantico acoplado a una microcavi- 
dad semiconductora" (Tesis de Dr. en Fisica, UdeA, en ejecucion). 

1.5.3 Direccion de proyectos de investigacion 

• Sistemas cuanticos de tres particulas, Universidad Nacional de Colombia - COLCIEN- 
CIAS, 1995-1996. 

• Sistemas cuanticos de pocas particulas en dos dimensiones, Universidad de Antioquia - 
COLCIENCIAS, 1997-1998. 

• Nucleos artificiales: sistemas de excitones en un punto cuantico, ICIMAF - CITMA, 
1999-2000. 

• Efectos magneto-opticos lineales y no lineales en puntos cuanticos semiconductores, ICI- 
MAF - CITMA, 2001-2003. 
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Tabla 1.1: Citas a los trabajos presentados de acuerdo con el ISI 



1.5.4 Reconocimientos 

• Premio Nacional de la Academia de Ciencias de Cuba en 1993 a la coleccion de trabajos 
denominada "Mecanica cuantica de pocos cuerpos y el metodo 1/D". 

• Premio Nacional de la Academia de Ciencias de Cuba en el 2001 a la coleccion de tra- 
bajos denominada "Espectro de energias, densidad de niveles, polarization del espin, 
propiedades de transporte y opticas de puntos cuanticos y trampas de atomos" . 

• Participation en el Premio Nacional de la Academia de Ciencias de Cuba en el 2002 a 
la coleccion de trabajos denominada "Dispersion Raman e hiper-Raman en sistemas de 
puntos cuanticos" , cuyo autor principal es E. Menendez Proupin. 

• Election de A. Gonzalez como Miembro Asociado del ICTP en el sector de Materia 
Condensada en el periodo de 1996 al 2003. 



1.5.5 Citas a los trabajos presentados 

Como un indice adicional del impacto de los trabajos realizados presentamos las citas a los 
mismos. Hemos separado las autocitas de las citas externas, es decir de las realizadas por autores 
sin ninguna relation con nosotros. Presentamos datos del Institute of Scientific Information 
(ISI). La information esta actualizada hasta Julio del 2005. 



Capitulo 2 



Descripcion tecnica abreviada 



2.1 El desarrollo en potencias de 1/D y su aplicacion en 
hadrones, conglomerados atomicos y puntos cuanticos 
semiconductores 

Como se menciono en la Introduction, el desarrollo en potencias de 1/D, donde D es la di- 
mension espacial o el momento angular, fue el metodo utilizado en el trabajo de candidatura 
del autor [UH] para obtener relaciones universales (independientes del potencial de interaccion) 
entre las magnitudes fisicas correspondientes a los sistemas de dos y tres particulas. De forma 
que los trabajos relacionados con la presente section El HI M UH1 UH1 EES] constituyen una 
extension directa del trabajo de candidatura. De ellos hemos seleccionado los articulos |2j, jl] 
y [12] j los cuales contienen algunos de los resultados mas interesantes y representativos en esta 
direction. 

La idea del metodo 1/D es sencilla y consiste en utilizar la dimension (o el momento angular) 
como parametro para organizar una solution en serie de la ecuacion de Schrodinger. El cuadro 
cualitativo que surge en la aproximacion principal del metodo es muy interesante. 

Consideremos la ecuacion de Schrodinger para dos particulas identicas con interaccion central 
en el sistema del centro de masa. La ecuacion radial en D dimensiones se escribe asi 



h 2 ( d 2 D-ld l(l + D-2)\ , 

lf + 1 5 " +gV(r)}^ = E^, (2.1) 



2m \dr 2 r dr r 2 

donde / es el mimero cuantico del momento angular y m = mo/2 es la masa reducida del par. 
Haciendo el cambio de variables ip = C/ r 2 j eliminamos el termino lineal en la primera 
derivada, lo cual da lugar al siguiente hamiltoniano: 

h 2 ( d 2 A(A + 1)\ , . 
h 2 = -—[ — - 1 2 ; + gV(r), (2.2) 

siendo A = I + {D — 3)/2. Tanto / — > oo como D — » oo son limites en los que la energia cinetica 
queda dominada por el termino centrifuge Tomando I = 0, por ejemplo, y D grandes resulta 
que el termino principal de h viene dado por la energia potencial efectiva: 
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U 2 (r) = ^ + gV(r). (2.3) 

Minimizando U 2 (r) hallamos un valor de r (una orbita), ro, y una aproximacion para la energia, 
^2(^0) • Correcciones a esta energia pueden hallarse escribiendo r = r + 5r y considerando 
oscilaciones (armonicas y anarmonicas) alrededor de r . Resulta: 

E 2 {g) = U 2 (r ) + (n + l/2)hw - + • • ■ , (2.4) 



donde 00(g) = ^/K 2 / (m 2 ^) + Agr^V'^r^/m/D es la frecuencia de pequenas oscilaciones. La 
derivada del potencial se toma respecto a r 2 . 

Expresiones analogas pueden hallarse para tres y, en general, para N particulas [HZj • En el caso 
de tres particulas identicas con interacciones por pares en estados con momento angular total 
igual a cero, resulta: 



ft 2 f 1 d 2 Id 2 ( 1 1 \ d 

h = -—<—^ + -^ + 



h 2 D 2 { 1 1 

+ 7r-^(^ + ^)+g(V 12 + V 23 + V 31 ) + ---, (2.5) 



2 \mdr 2 pdp 2 \mr 2 pp 2 J d6 2 
2 D 2 ( 1 J_ 
8 sin 2 9 \mr 2 pp 2 

donde r = r 2 — r\, p = r 3 — (r*i + f 2 )/2, cos 9 = r ■ pj (rp) y las masas reducidas son m = mo/2, 
p = 2m /3. 

En el limite D — > 00, la energia potencial efectiva que sobrevive en el hamiltoniano, U 3 , es min- 
imizada en una configuracion del tipo triangulo equilatero (cos 9 = 0, p 2 = 3r 2 /4), resultando: 

U 3 (r, g) = ^ + 3gV(r) = 2U 2 (r, 3g/2), (2.6) 

por lo que |o"7] : 

E 3 (g) 



2E 2 (3g/2) 



= 1. (2.7) 



Esta relacion se satisface en muchos casos con un error inferior al 10 % incluso en el limite fisico 
D = 3. Notar que es una relacion universal, es decir independiente del potencial de interaction, 
V. 

Correcciones a la Ec. ()2.6|) se calculan relajando la configuracion rigida del limite D —>■ 00, 
resultando: 



E 3 (n+, n_,n e ) = 2U 2 (r , 3g/2) + (n + + l/2)tvuj+ + (n_ + l/2)^_ 



3ft 2 



+ (^ + 1/2)^- (2i 
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donde ahora existen tres frecuencias de oscilaciones, correspondientes a un modo simetrico 
( "respiratorio" ) , u + , a un modo asimetrico, w_ y a un modo de deformation, ug. Interesante 
el hecho de que: 



!ko + (g) = hu(3g/2), (2.9) 

lo cual relaciona la energia de excitation del modo simetrico con la energia de excitation del 
sistema de dos particulas. Esta es tambien una relation universal. 



2.1.1 Razon entre los radios cuadraticos medios de bariones y mesones 

Comparando las secciones elastica e inelastica en la dispersion hadron-proton, Povh y Hufner 
[HH1 inH HOj estimaron los siguientes valores para los radios cuadraticos medios de mesones y 
bariones (en fm 2 ): 



# 2 =0.41, # 2 W = 0.52, #^ = 0.35, #5 = 0.21, (2.10) 
R 2 ^ = 0.04, #3 = 0.67, #i )E = 0.58, #2=0.50. (2.11) 

En el trabajo [2] se hizo un analisis de estos valores a partir del metodo 1/D. En particular 
se analizo la relation #bar y on/-^meson- El trabajo contimia la direction iniciada en jH] donde se 
obtuvieron relaciones entre las masas de baryones y mesones. 

En la aproximacion principal del metodo 1/D, la distancia cuadratica media entre las particulas 
de un meson, compuesto de un quark q y el correspondiente antiquark q, se determina de: 

(r 2 (g)) q , = r 2 + O(l/D), (2.12) 
donde r es la distancia que minimiza el potential efectivo: 

fi r = m q /2 es la masa reducida del par y hemos separado explicit amente la constante de 
acoplamiento g. Al final haremos g — 1. 

En el baryon, compuesto de tres quarks, el potential efectivo en el limite D — » oo depende de 
dos coordenadas de Jacobi: r y p. En estados con momento angular total igual a cero, solo las 
combinaciones escalares r 2 , p 2 y 7 = r • p/{rp) entran el potential. Cuando los tres quarks son 
identicos, el minimo se alcanza en una configuration del tipo triangulo equilatero, resultando: 

U qqq (g) = 2U qq (3g/4), (2.14) 
por lo que para las distancias obtenemos: 

(r 2 (9)) qqq =(r 2 (3g/A)) qq + 0(1/D). (2.15) 
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Las ecuaciones (|2.6j) y (|2.14|) se diferencian en el argumento. En la Ec. (|2.14|) el potencial de 
un par esta evaluado en 3g/4, mientras que en la (|2.fij) esta evaluado en 3g/2. La razon es que 
el potencial de interaccion entre dos quarks en el baryon se relaciona con el potencial entre un 
quark y un antiquark en el meson asi: 

Vqq\baryon ~ ~^^qq\mesori' (2.16) 

Los radios cuadraticos medios, definidos como valores medios de las distancias al centro de 
masa al cuadrado, satisfaran, por tanto: 

R qq -{3g/4:) 3 

Para comparar directamente los radios del baryon y el meson (es decir R^(g)) es necesario una 
hipotesis sobre el potencial de interaccion, V. Utilizaremos el potencial de Martin [72], para el 
cual V(r) ~ r@, con (3 = 0.1. Obtenemos fmalmente: 



D2 //L n (/9+4)/08+2) 



Rqq V • > 



+ 0(1/D). (2.18) 



1.40 q 




0,90 'T* rrr * t " T " 1 ' ' ' ' M' ' m 1 1 1 t i \ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i i n 
0,0 1.0 2.0 3.0 4,0 



Figura 2.1: Comparacion entre los resultados del metodo 1/D (linea solida) con los valores 
experimentales para el cociente -Rbaryon/ -Rmeson- El valor predicho por el modelo naif de partones 
se da como referenda. 



En el caso de los quarks mas ligeros (u, d y s) es necesario incluir correcciones al termino 
principal que provienen de la relatividad del movimiento. Expresiones analogas a ()2.18j) pueden 
obtenerse para mesones compuestos de quarks distintos (qQ y Qq) y baryones del tipo Qqq. Los 
resultados para el cociente -Rbaryon/ -Rmeson se muestran en la Fig. l2.1l como funcion del parametro 
£ = mQ/m q . Como referenda, se muestra el valor predicho por el modelo "naive" de partones, 



CAPITULO 2. DESCRIPCION TECNICA ABREVIADA 



22 



a/3/2. Aunque los valores medidos de -Rbar y on/-Rmeson no son conclusivos, la curva obtenida en 
la aproximacion principal del metodo 1/D reproduce cualitativamente la dependencia con £. 

2.1.2 Energias de ligadura y de excitacion en clusters atomicos con 
interaccion de Lennard-Jones 

En el trabajo [I] se estudiaron conglomerados de atomos de gases inertes (Ne, Ar, Kr, Xe) a 
partir del metodo 1/D. Para estos atomos, el potencial de interaccion puede aproximarse con 
un potencial central. Nosotros utilizamos la parametrizacion de Lennard-Jones 12-6 dada en el 
articulo [TSj- Los conglomerados atomicos han sido estudiados extensivamente en los ultimos 
anos. Los sistemas mas pequenos (B3, N4, etc.) han sido analizados por metodos ab-initio 
(ver, por ej., jHHZSj), mientras que en los mas grandes (cientos de atomos) se ha utilizado una 
combination de la dinamica molecular con el metodo de los funcionales de la densidad [To] . 

La propuesta nuestra para estimar la energia del estado base del cluster parte de la relacion 
universal (independiente del potencial de interaccion) obtenida en la aproximacion principal 
del metodo 1/D entre las energias del complejo de n particulas y de dos particulas [HHl EJJ : 
E n (g) « (n — l)E 2 (ng/2), donde hemos escrito el potencial de interaccion entre dos particulas 
como gV, separando explicit amente la constante de acoplamiento g. Al final escribiremos 
g = 1. Correcciones a esta relacion pueden ser explicitamente calculadas con el metodo 1/D, 
result ando: 



R = Eni9) -1 I 1 1 f H '\- 

(n-l)E 2 {ng/2) D (n - l)e 2 (ng/2) \ 2 

n—1 n(n — 3) (n— l)(n — 2)1 M , . „ 9 . 
+ -^ + -4-^- ^ 2 ^ >- \+0(l/D>). (2.19) 

En esta formula, e 2 es el termino principal en el limite D — ► 00 de la energia del dimer (dos 
atomos), u , uj\ y u 2 son las frecuencias de pequenas oscilaciones del cluster de n atomos 
alrededor de la estructura rigida que se obtiene en el limite D — >• 00, mientras que \i\ es la masa 
reducida del dimer y p x la distancia media entre dos atomos del cluster. 

La formula (J2.19|) debe complementarse con un factor geometrico ya que la interaccion entre 
atomos es de corto alcance . En D — » 00 cada atomo tiene como vecinos a todos los demas, 
pero en D = 3 no. La ecuacion final es: 

/ min V . . . V(fi — r,-) \ 
E n {g) = . VM , R(n - l)E 2 (ng/2). (2.20) 

\n(n — 1) rmnV(rj/2 / 

El factor entre parentesis tiene el sentido de cociente entre el numero efectivo de interacciones 
por pares en D = 3 y D = 00. Con esta expresion evaluamos la energia de clusters de n atomos 
(donde n va desde 3 hasta 13) y las comparamos con los resultados por Monte Carlo de la 
Ref . [73] . Los resultados para el Ar se muestran en el Cuadro 12.11 En el peor de los casos, la 
diferencia entre ambos es del orden del 1 %. 

Por otro lado, comparando las frecuencias de pequenas oscilaciones en el trimer (tres atomos) 
y el dimer podemos estimar la energia de excitacion del primero en terminos de la del segundo, 
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n 


MC 


l/D 


3 


-2.553 


-2.558 


4 


-5.113 


-5.131 


5 


-7.791 


-7.804 


6 


-10.88 


-10.91 


7 


-14.18 


-14.19 


13 


-38.65 


-38.35 



Tabla 2.1: Comparacion entre los resultados por Monte Carlo y el metodo l/D para la energia 
de clusters de n atomos de Ar. 
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Figura 2.2: Comparacion entre los resultados del metodo l/D (linea solida) con los valores 
numericos de la Ref. para la energia de excitacion de los trimers. 




result ando: 



AE 3 (g) = AE 2 (3g/2) {1 + 0(1/D)} . (2.21) 

En la Fig. 12.21 mostramos la comparacion de nuestro estimado con los calculos numericos del 
articulo [7^]. Vemos que el error es del orden del 30% en el Ne 3 (el cluster mas ligero), donde 
los efectos cuanticos deben ser mayores, pero solo del 3% en el Xe 3 . 



A 

2.1.3 Espectro de energias de pocos electrones en un punto cuantico 

En el trabajo [THj se estudiaron los niveles de energia de dos y tres electrones en un punto 
cuantico cuasi-bidimensional. Sobre los electrones actiia un campo magnetico ortogonal al 
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piano del movimiento. El espiritu del articulo es muy similar al de la Ref. [77], donde se 
hace un analisis cuasi-clasico del mismo problema. La diferencia radica en el parametro del 
desarrollo en serie, que es 1/\J\ en nuestro caso (J es el momento angular total) y h en la Ref. 

E3 

El punto cuantico se modela con un potencial parabolico en el piano. En este caso el movimiento 
del centro de masa se puede separar y para el problema de N electrones el hamiltoniano escalado 
que describe el movimiento relativo es: 



h 




(2.22) 



donde /4 son los correspondientes vectores de Jacobi en el piano, 8k es el angulo entre p k y Pk+i, 
(3 = (E cou i/ {hVt)) l l & , donde E cou \ es la energia caracteristica de Coulomb y Q = a/^o + uj 2 /4, 
siendo fruo la energia del confinamiento y u c la frecuencia ciclotronica en el campo B. La 
constante (3 se ha "renormalizado" , /3 3 = /3 3 /|J| 3//2 , de forma que en el limite formal \J\ — > oo 
la constante (3 se tomara fija y Coulomb se tratara no perturbativamente. En la Ec. ()2.22|) las 
longitudes se han escalado de acuerdo a p 2 — > \J\p 2 . 

Puede verse entonces que solo la energia potencial efectiva (centrifuga + confinamiento + 
Coulomb) sobrevive en la Ec. ()2.22|) en el limite \J\ — > oo. Las primeras correcciones provienen 
del movimiento armonico de pequena amplitud alrededor de la configuracion de equilibrio (el 
minimo del potencial efectivo). Correcciones mas altas provienen de considerar la anarmonici- 
dad. 

En el trabajo [13] se calcularon los coeficientes del desarrollo hasta el para el caso de dos 
electrones y hasta el para tres. La serie para e (es decir, el autovalor de h) se escribe: 



Las comparaciones de nuestros estimados con calculos numericos exactos para dos electrones 
nos permiten concluir que el error al estimar la energia no rebasa el 0.3%. 

En el caso de tres electrones podemos comparar la serie en 1/\J\ con los aproximantes de Pade 
obtenidos en JT]. Un ejemplo de esta comparacion se muestra en la Fig. 12.31 Vemos que tanto 
el Pade como la serie en 1/\J\ son asintoticamente exactos en los limites (3 — > y (3 —>■ oo. La 
diferencia entre ambos estimados no rebasa el 0.4% en todo el rango de variacion de (3. 



CAPITULO 2. DESCRIPCION TECNICA ABREVIADA 



25 



\PtA-\JU\ 
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Figura 2.3: Diferencia entre la serie en 1/\J\ y el aproximante de Pade para la energia del 
estado de tres electrones con numeros cuanticos J = —3 j S = 3/2. 

2.2 Cuantizacion semiclasica y supersimetria en Meca- 
nica Cuantica. Vibraciones moleculares 

Como se menciono en la section a partir de la supersimetria en el articulo pQ se construyeron 
potenciales de interaction entre N particulas en una dimension para los cuales se conoce ex- 
actamente la energia y la funcion de onda del estado base, mientras que la supersimetria en 
conjuncion con la cuantizacion semiclasica y las transformaciones canonicas se utilizo en los 
trabajos 13 IE] para obtener los estados excitados de moleculas lineales y los estados de 
tres particulas con momento angular total igual a cero. 

En la presente section nos concentraremos en la mejora de la description semiclasica del espectro 
de excitaciones en sistemas de tres particulas, lo cual nos parece el resultado mas interesante 
de los trabajos mencionados. En la Ref. [Hj se argumento que reemplazar el potential de 
interaction, V, por un nuevo potential, U, construido a partir de la funcion de onda del estado 
base, debia conducir a una mejor description semiclasica de los estados excitados en un sistema 
con varios grados de libertad. La relation entre V y U es la siguiente. Sea 

H = lL^r + v > ( 2 - 24 ) 

el hamiltoniano y \I/o) E la funcion de onda y energia del estado base, respectivamente. V 
puede ser escrito en terminos de como: 



V-E = J2^{(d t Wy-hd*W}, 



(2.25) 
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donde = ex P —W/h. Entonces, el nuevo potencial U es definido de la siguiente forma: 



U = J2^-(d l W) 2 . (2.26) 

Listemos algunas de las propiedades de U. 

- U es definido positivo. U = en el minimo, que corresponde a la configuration de maxima 
probabilidad, es decir, el maximo de |^o| 2 - Notemos que E' = es la energia exacta del 
estado base en el potencial V — E . 

- En el caso de un grado de libertad, la regla de Bohr-Sommerfeld modificada: 



^ j> ^2fi(AE - Ujdx = nh, (2.27) 

resulta exacta en una clase de potenciales que incluye el potencial de Coulomb Estos 
potenciales tienen la propiedad de que las parejas supersimetricas 

H± - E = — {p 2 + (dW) 2 ± hd 2 W\ , (2.28) 
2/x k J 

son invariantes de forma, es decir tienen la misma dependencia funcional y difieren solo 
en los parametros de que dependen. Existen indicaciones j3] de que una clase similar 
existe en sistemas con varios grados de libertad. 

- Formalmente, puede pensarse que U difiere de V—E en un termino que es de orden h. Por 
eso en el limite de numeros cuanticos muy grandes debe ser lo mismo hacer la semiclasica 
con U o con V — Eq. Procediendo formalmente es posible argumentar que la regla de 
Bohr - Sommerfeld se transforma en la Ec. ([2.27)1 . Y esta ultima tiene la propiedad de ser 
exacta tambien en el limite n = ya que el estado base tiene exactamente E' Q = 0. Por 
eso los niveles de energia obtenidos a partir de la Ec. (|2.27|) son asintoticamente exactos 
cuando n -> oo y exactos por construccion cuando n — 0. De aqui que la semiclasica en 
U sea mejor. 



En el caso de varios grados de libertad, las reglas de cuantizacion a utilizar son las de Einstein 
- Keller - Maslov jTHlIZi: 

J2Pi dx i = K + m a /A)h, (2.29) 

' a i 

donde V — E se sustituye por U, los indices de Maslov se ponen a cero m a = y el estado 
base corresponde a todos los n a = 0. 

Debido a la forma funcional de U, con un minimo absoluto U — 0, las energias de excitacion mas 
bajas pueden calcularse aproximadamente a partir de las frecuencias de pequenas oscilaciones 
alrededor del minimo 0. Esta descripcion corresponde con el cuadro semiclasico de los estados 
excitados como oscilaciones por encima del estado base. 
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Una description mas completa requiere tomar en cuenta la anarmonicidad de U cerca del 
minimo. Para ello utilizamos la reduction a formas normales de Birkhoff y Gustavson (BG) en 
la cual la anarmonicidad es tratada no perturbativamente por medio de sucesivas transforma- 
ciones canonicas [7^]. Esquematicamente, el algoritmo se puede formular asi. Desarrollando 
el potential alrededor de la configuration de minimo, x[°\ obtenemos una expresion para el 
hamiltoniano clasico: 

H(p, q) = V{x<®) + H 2 (p, q) + H 3 (p, q) + • • • , (2.30) 

donde los H a son polinomios homogeneos de grado a en las variables pi y qi = x. t — xf^ . Ahora 
hacemos transfer maciones canonicas sucesivas (la forma explicit a de las cuales puede verse en 
|79| |H]) para reducir H 2 , H 3 , etc a la forma normal (la forma normal de H 3 y de todos los 
polinomios impares es cero). De forma que el hamiltoniano resulta: 

T(P, Q) = V{xV>) + \Y^< P ! + Qi) + T 2(P 2 + Q 2 ) + UP 2 + Q 2 ) + ■ ■ ■ , (2.31) 

i 

donde los T a son polinomios homogeneos en las variables P 2 +Q 2 . Truncando la expansion (|2.31|) 
obtenemos un hamiltoniano que describe aproximadamente el problema y que es exactamente 
integrable. La receta de cuantizacion es simplemente: 

(P 2 + Q 2 )/2-+(n t + l/2)h, (2.32) 
o, en el caso de que el punto de partida es el hamiltoniano modificado, con el potential U, 

(P^ + Ql)/2^ ni h. (2.33) 




Figura 2.4: Energia de excitation del estado con numeros cuanticos (n 1 ,ra 2 ) = (1,1) en el 
modelo de molecula triatomica lineal del articulo [7j. 
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En la Fig. 12.41 mostramos un ejemplo de como trabaja la semiclasica modificada en un modelo 
de molecula triatomica en una dimension [2j . Separando el centro de masa restan dos grados de 
libertad, es decir dos numeros cuanticos ni, n 2 . El parametro A cuantifica la anarmonicidad de 
la molecula. La figura muestra la energia de excitacion del estado con numeros (1, 1). La curva 
denominada BG corresponde a la cuantizacion a traves de BG del movimiento en el potencial V, 
mientras que BG* corresponde al movimiento en U. Los resultados se comparan con estimados 
numericos utilizando una base de osciladores. Los resultados de BG* son significativamente 
superiores cuando la anarmonicidad es grande. 




aeon 1 1 1 1 ri ti 

0,1 0.2 OA 0.6 0-6 

lambda 

Figura 2.5: Comparacion entre BG* (linea solida) y los estimados numericos (cuadrados) para 
los modos de oscilacion del tipo "respiratorios" en el modelo de molecula triatomica del articulo 
0. 



Un modelo de molecula triatomica en tres dimensiones es estudiado en 0. Se consideran 
solo estados en los que el momento angular total es cero. Tambien existe un parametro A que 
caracteriza la anarmonicidad del potencial. En la Fig. I2.5l se muestran las energias de excitacion 
de los modos "respiratorios" (es decir donde la molecula se dilata y contrae a lo largo de la 
direction radial) como funcion de A y se comparan con estimados numericos. La concordancia 
tambien es muy buena, siendo el error maximo del orden de 1%. 

2.3 El metodo de los aproximantes dobles de Pade en 
sistemas con decenas de electrones o de atomos con- 
finados 

El metodo de los aproximantes de Pade de dos puntos fue utilizado por primera vez en los 
trabajos fH^JEJCElEl P ara obtener aproximaciones analiticas a la energia del estado base 
y los primeros estados excitados de electrones en puntos cuanticos y atomos con espin entero 
(bosones) confinados en trampas. En esta seccion resenaremos los trabajos [TI], ^3] y |16j . 
considerados los mas importantes. 
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Un aproximante de Pade de dos puntos es una funcion racional de la variable (3 que interpola 
entre dos series conocidas alrededor de los puntos (3 — (3q y /3 = (3\. En el caso nuestro, (5q — 
y la serie alrededor de j3 = no es mas que la teoria de perturbaciones hasta determinado 
orden. Por otro lado, f3\ = oo y la serie de Laurent alrededor de (3x tambien se calculara hasta 
determinado orden. En general, obtendremos los desarrollos siguientes: 



e|^ = $> fe /3 fc + 0O3 s+1 ), (2.34) 

k=0 

e\^oo = (3 2 lj2a k /p k + 0(l/p t+1 )\. (2.35) 
U=o J 



El aproximante doble de Pade se escribira: 



Como hemos asumido e^-^ ~ /3 2 , entonces L = K + 2. Sin perder generalidad, podemos hacer 
go = 1- La igualdad entre los mimeros de coeficientes p, q por un lado y b, a por otro resulta en 
que: 



s + t = 2K+l. (2.37) 

Igualando las ecuaciones ()2.34p2.35j) con ()2.3(ij) obtenemos un sistema de ecuaciones lineales 
para los coeficientes Pk J Qk- Los aproximantes que mejores propiedades de convergencia tienen 
son aquellos en que s ~ t. 



2.3.1 Energia del estado base de puntos cuanticos con decenas o 
cientos de electrones 

____ / ^ 

En el trabajo [IT] se elaboraron aproximaciones analiticas a la energia de sistemas de hasta 
5 electrones en un punto cuantico cuasi-bidimensional y en presencia de un campo magnetico 
externo orientado perpendicularmente al piano del punto. La idea original se pernio en el 
trabajo [53] donde se calcularon los niveles del exciton y se basa en el hecho de que en el 
limite de altas densidades (confinamiento fuerte) la energia se puede calcular por teoria de 
perturbaciones, mientras que en el limite opuesto de densidades muy bajas (confinamiento 
debil) el sistema forma una especie de solido (molecula) de Wigner cuya energia de Madelung 
y frecuencias de pequenas oscilaciones tambien pueden ser calculadas. La interpolacion entre 
ambos regimenes se hace utilizando un aproximante de Pade de dos puntos. La idea habia sido 
empleada en los 70 para estimar la energia del gas de electrones a densidades intermedias [80j. 
Posteriormente, el calculo se extendio a puntos cuanticos con decenas o cientos de electrones 

El punto de partida es el hamiltoniano del sistema escrito en variables adimensionales: 
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TT 1 N 1 

donde u; es la frecuencia del punto y el parametro j3 se define en terminos de las energias 
caracteristicas de Coulomb y de oscilador, /3 3 = ^E cou i/(fkoo). El limite (3 — > corresponde a 
electrones libres, mientras que /3 — ■> oo corresponde a la molecula de Wigner. Las series en el 
entorno de estos puntos se escriben: 



e|^o = b + b 3 p 3 + ..., (2.39) 
= (3 2 {a + a 2 /(3 2 + ...}. (2.40) 

Los coeficientes b a y a a fueron calculados para puntos con un maximo de 210 electrones. b 
es la energia de N electrones libres en un potencial armonico, mientras que 6 3 coresponde al 
calculo del valor medio de la energia de interaccion Coulombiana en primer orden de teoria de 
perturbaciones. Por otro lado, escalando las distancias segiin r — > j3p, podemos reescribir h asi: 



h 



Por lo que, en el limite (3 — > oo, la energia se determina de minimizar la energia potencial 
efectiva: 




(2.42) 



(una especie de energia de Madelung para un sistema finite) , mientras que las correcciones se cal- 
culan a partir de las oscilaciones de la estructura que minimiza U. En el caso N < 5 uno puede 
adivinar la geometria del minimo. Por ejemplo, en el caso iV = 3 uno busca el minimo como un 
triangulo equilatero. Pero para iV mas grandes la minimizacion debe hacerse numericamente. 
Nosotros utilizamos una biisqueda aleatoria que despues es optimizada con metodos de gradi- 
ente. Los modos normales de oscilacion tambien fueron obtenidos numericamente. 

El comportamiento asintotico cuando el numero de electrones es grande, iV > 100, es el sigu- 
iente: b ss |A^ 3//2 , 6 3 ss 0.7 iV 7 / 4 , a w 1.0 iV 5 / 3 y a 2 w 0.6 N 5 ^. 

Los aproximantes dobles de Pade a utilizar son: 



PsM = b + a (3 2 l\ 
PaM = b + 



1 



l + (b 3 /a )f3 + (a /(b -a 2 ))f3* 



(2.43) 



+ a p z ->l- 



l + qiP + g 2 /3 2 + q 3 (3 3 



l + q l f3 + q 2 ^ + q^\ ' 



(2.44) 
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Figura 2.6: Los limites (3 —>■ y j3 —>■ oo (curvas discontinuas) y el aproximante P 4j3 (curva 
continua) para el estado base del punto con 42 electrones. 



Los coeficientes q a que aparecen en la Ec. (|2.44jl son determinados de la condition que los 
desarrollos en/3^0y/3^oo sean reproducidos. 

En la Fig. 12.61 mostramos el comportamiento del Pade en el punto con 42 electrones. El 
estado base tiene momento angular total y espin total L = y S = 0, respectivamente. El 
aproximante interpola suavemente entre los dos regimenes. El error estimado no es mayor que 
el 2.5 %. Notese que existe un intervalo de "densidades" (alrededor de j3 = 0.9) donde los dos 
desarrollos asintoticos mantienen su validez. 



2.3.2 Aproximaciones analiticas a la solucion de la ecuacion de Gross- 
Pitaevskii para bosones en trampas atomicas 

En el trabajo se utilizaron los aproximantes dobles de Pade para obtener estimados analiticos 
a la energia y el potencial quimico de atomos con espin entero (bosones) en trampas atomicas 
a temperaturas practicamente iguales a K. Estos sistemas son usualmente descritos por la 
ecuacion de Gross-Pitaevskii (GP) [HE] la cual, en la variante mas simple de potencial externo, 
se escribe de forma adimensional asi: 



-A + -r 2 + g\ij\ 2 -n}ip = 



(2.45) 



CAPITULO 2. DESCRIPCION TECNICA ABREVIADA 



32 



3.5 -r 



3-- 



25 



2-. 



1.5 



0.5 -■ 




□ Ret 12 

— P«(R) 



4- 



1 

2.5 3 

R 



0.5 



1.5 



Figura 2.7: Comparacion entre el aproximante Pg^ y los resultados numericos del articulo 
para la energia por particula en dos dimensiones. 



donde if) es la funcion de onda del condensado de bosones, \i es el potencial quimico y g la 
constante de interaction. El limite g — > corresponde a bosones libres y las correcciones al 
mismo pueden calcularse por teoria de perturbaciones, mientras que en g — > oo la denominada 
teoria de Thomas-Fermi es aplicable [SHj y las correcciones a la misma son debidas a efectos en 
la frontera del condensado. En tres dimensiones, estas correcciones fueron calculadas por Fetter 
y Feder |81j . Calculos analogos se presentaron en la Ref. JIB] para el problema bidimensional. 
Notese que la Ec. (|2.45jl tiene la forma de una ecuacion de Schrodinger no lineal donde \i juega 
el papel del autovalor. 

En terminos de if), el potencial quimico y la energia del sistema de atomos se escriben: 



E 



d''r{r(-^ + \r 2 )iP + g\iP\ 4 



9 



(2.46) 
(2.47) 



Para estas magnitudes podemos obtener series en los limites g — > y g — >• oo: 



f(R)\ 



b + hR b + 0{R 10 ) 



(2.48) 
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E 


u 

1 


bo 


1 


1 




1/16 


1/8 


^8 




-0.001124 


-0.003371 


Q>o 


1/3 


1/2 


04 


4/3 


2/3 


Q2 


0.646107 


0.751621 


Qi 


0.292455 


0.424002 




0.073825 


0.142145 


P4 


0.262488 


0.379052 



Tabla 2.2: Coeficientes de los desarrollos asintoticos para la energia y el potencial quimico de 
bosones en dos dimensiones. 



f{R)\ R ^ = R 2 {a + ^H AR ) + °( 1 / R5 )}^ ( 2 - 49 ) 

donde R = (15(?/(47r)) 1 / 5 tiene interpretation de radio del condensado en el limite de g grandes. 
A partir de las mismas podemos construir los aproximantes de Pade, por ej. 



Q3 = h/a , q% = 3g 3 /(2a ), 

<?4 = b^/ao, q 2 = 3b 5 q 1 /(2a 2 ), (2.50) 

los cuales brindan aproximaciones analiticas para estas magnitudes a cualquier R. El hecho de 
que el termino que acompana a CI4 contenga la dependencia ln(_R), que no es polinomica, puede 
resolverse construyendo el aproximante para la magnitud (l/R)d/dR (R 2 f(R)). 

En la Fig. l2.7l mostramos como se compara el aproximante P§^(R) con la solucion numerica de 
la ecuacion de GP obtenida en el trabajo [H2| para la energia de bosones en dos dimensiones. 
Ahora R = (Ag/n) 1 ^ y el aproximante se define segun: 

p h f*A 5 l + g 2 x 2 + P 4 x 4 f 



R 2 Jo 1 + Q2X 2 + Q4X 4 + Q 6 x£ 

donde los coeficientes que en el intervienen vienen dados en el CuadroEZO La diferencia entre 
el Pade y la solucion numerica no excede el 2 % en el intervalo de R considerado. 



2.4 La funcion BCS y las ecuaciones de Bethe-Goldstone 
para sistemas multiexcitonicos confinados 

En los trabajos [THl 123 EH3 EH] se hace uso extensivo de la funcion BCS para describir el 
apareamiento de fermiones, mientras que en el articulo [20] se utilizo la denominada ecuacion 
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de Bethe - Goldstone para estudiar sistemas de hasta 12 pares electron - hueco y sistemas no 
neutros (hasta 6 electrones y 2 huecos). Mas abajo describimos los resultados mas relevantes 
de estos articulos. 



2.4.1 La transicion al regimen BCS, la luminiscencia coherente y 
otras propiedades de sistemas multiexcitonicos 

En el tr abajo ^Hj se resalto e hizo uso de la analogfa entre los micleos atomicos y los sistemas 
de excitones en semiconductores. Un metodo de calculo de la estructura nuclear muy empleado 
para micleos pesados y que se basa en una funcion BCS [SHI fue utilizado para calcular la 
energfa de sistemas de hasta 90 pares de electrones y huecos, es decir hasta 180 partfculas, con 
el objetivo de investigar cuando, o sea para que densidades, el sistema se comporta como un 
gas de electrones y huecos y cuando se manifiesta como un sistema de excitones descrito por la 
funcion BCS. 

La funcion de onda variacional BCS, 



\BCS) = H (uj + Vjhtetm, (2.52) 

aparea huecos y electrones en estados con momento angular total igual a cero. El estado j del 
electron tiene momento invertido respecto del hueco. \uj\ 2 es la probabilidad de que no exista 
el par (j, j), mientras que \vj\ 2 es la probabilidad de que exista. El principio de Ritz 



E gs < (BCS\H\BCS), (2.53) 

nos da un estimado para la energia del estado base, E gs . Los coeficientes Uj y Vj se utilizan 
como parametros variacionales sujetos a la ligadura \uj\ 2 + \vj\ 2 = 1. La minimization conduce 
a la denominada ecuacion del gap 

A,- = y~7j, k\V\kJ) — Ak =, (2.54) 

donde V es el potencial de interaction entre pares, e^ F es la energfa del estado k renormalizada 
por la interaction con el resto de los estados y es el potencial qufmico que se determina de la 
ecuacion 

N= (BCS\J24 e i\ BCS )- ( 2 - 55 ) 

j 

Los coeficientes Vj se han escrito en terminos de los parametros de gap, Aj, asf: 
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Las ecuaciones (|2.54j) se resuelven iterativamente y muestran muy buenas propiedades de con- 
vergencia. 



1.5 
1.25 



-0.25 
-0.5 





O 6 particles 




□ 1 2 particles 




20 particles 




A 30 particles 




* 42 particles 


- 


%^ 4- 56 particles 




^» x 72 particles 




V 90 particles 




HF 










- 










% \ 




X ^ 




x ^ 




K 




o 








0.25 0.5 0.75 1 1.25 

P/N 1 ' 4 



Figura 2.8: Escalamiento de la energia en sistemas de hasta 90 pares electron - hueco. La 
curva continua representa el estimado de Hartree - Fock. 

Debido a que la funcion BCS no conserva el mimero de particulas se hace necesario un pro- 
cedimiento de proyeccion. En el trabajo se empleo una proyeccion aproximada propuesta por 
Lipkin y Nogami [SHI El] • Al resultar evidente que esta proyeccion era insuficiente se creo un 
metodo estocastico basado en el algoritmo de Metropolis que ofrecio muy buenos resultados, 
incluso en sistemas nucleares [TH] . 

Ademas del estimado de la energia de iV excitones (2 < < 90) para cualquier valor del 
confinamiento del punto cuantico, las conclusiones mas importantes del articulo JH] son las 
siguientes: a) que la energia muestra propiedades de escalamiento y b) que el apareamiento 
es significativo cuando (3 > 0.55 iV 1 / 4 . Traducido a densidades, esto significa que la densidad 
superficial de pares es mayor que 4/(7ia 2 B ), donde a# es el radio de Bohr efectivo del material. 
Estas propiedades son evidentes en la fig. 12 .81 El apareamiento aqui se visualiza como separacion 
de los valores de energia respecto del estimado de Hartree - Fock. 

Por otro lado, en el trabajo se analizaron los efectos combinados de la energia de con- 
finamiento del punto cuantico, la energia Zeeman, la interaction de Coulomb y de campos 
magneticos intensos sobre la polarization de los espines electronicos y la luminiscencia coher- 
ente en sistemas de hasta 40 pares e - h. 

El esquema teorico ya habia sido utilizado en el articulo [12], es decir una funcion de onda BCS 
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con proyeccion aproximada a lo Lipkin - Nogami (LN). La bondad de este metodo en el caso de 
campos magneticos grandes consiste en que la funcion BCS reproduce el estado basico exacto 
- un condensado de excitones [EH]- en el lmiite de campos infinites. 




30 35 
Magnetic Field [ T ; 



Figura 2.9: Position e intensidad de la linea de luminiscencia como funcion del campo 
magnetico y el numero de excitones en el punto. Ex es la energia del exciten. 



Mediciones de la luminiscencia en pozos cuanticos como funcion de la intensidad y polarizacion 
de la luz que crea los excitones habfan sido reportadas en 

Nuestros calculos arrojan solo un 10 % de polarizacion neta de los espines electronicos, aiin 
a campos tan altos como 20 Teslas, debido a la fuerte interaccion entre electrones y huecos. 
Ademas, se predice un corrimiento hacia el azul del pico de luminiscencia coherente con el 
aumento del numero de pares e - h presentes en el punto (o, lo que es equivalente, aumento de 
la intensidad del laser). Un ejemplo de estos resultados se muestra en la Fig. 12.91 
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2.4.2 Dinamica de un punto cuantico acoplado a una microcavidad 
optica 

En el trabajo [HE] se estudia la dinamica de un punto cuantico acoplado resonantemente a una 
microcavidad optica. La motivation es un trabajo teorico previo jHZ| que a su vez se inspira 
en resultados experimentales recientes [SHI EH E2] • En |HZ] se describe el campo fotonico en 
terminos de un espacio de Fock que contiene hasta 100 fotones. Considerando un solo estado 
del exciton en el punto, la dinamica en esta base de funciones se describe por un sistema de 
300 ecuaciones lineales acopladas. En el trabajo nuestro, por el contrario, se utiliza para los 
fotones un campo coherente (descrito solamente por una amplitud y una fase) por lo que, aun 
siendo una description simplificada, permite incluir mayor cantidad de estados excitonicos en 
el punto y asi considerar otros efectos. 




Nuestro esquema de trabajo se basa en considerar la evolucion temporal de magnitudes como: 
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Pn = (4e n ) = (h{h n ), 
N ph = \a\ 2 = (a f a), 



(2.57) 
(2.58) 
(2.59) 



p n describe la ocupacion del nivel electronico n en el punto (o del nivel correspondiente n del 
hueco), a es el campo fotonico y K nn es la funcion de apareamiento o polarizacion. p n y K nn 
provienen de una descripcion del subsistema electronico por medio de una funcion BCS. 

Las ecuaciones dinamicas se hallan de calcular el conmutador con el hamiltoniano. Por ejemplo, 
para p n : 




A (ps ) 



Figura 2.11: El numero de fotones como funcion de A para el punto con N states = 2 y 
parametros P = 1 ps -1 , k = 0.1 ps~ x . En el panel superior, las interacciones de Coulomb no 
han sido incluidas ((3 = 0). 



(2.60) 
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Las mismas se complement an con terminos que describen fenomenologicamente las perdidas a 
traves de los espejos de la microcavidad, el decaimiento espontaneo de los estados excitonicos 
en el punto y el bombeo incoherente del punto cuantico. Por ejemplo, para p n se obtiene: 

^-K n ^2(n, j \1/ r\j,n)Kj sin (fa -fa) - K n sin(0-0„) 

1Pn + P{l-p n ). (2.61) 

En esta ecuacion, (3 es la constante de la interaccion de Coulomb, K n y fa son el modulo y la 
fase de K nn , mientras que s y <fi son el modulo y la fase de a. Las constantes 7 y P parametrizan 
respectivamente el decaimiento espontaneo y el bombeo incoherente del punto. Las perdidas 
de fotones en la microcavidad se caracterizan por la constante k. 

En la Fig. I2.1UI mostramos la evolution temporal de las magnitudes p\, K\, N ph y sin(0 — fa) 
para el caso en que en el punto hay un solo estado del par. Estos resultados se corresponden 
cualitativamente con los del trabajo [87j . 

Por el contrario, en la Fig. 12.111 mostramos el caso en que en el punto hay dos estados ex- 
citonicos. La magnitud ploteada es el numero de fotones, N p h, para tiempos largos como 
funcion de A, la cual tiene interpretation de diferencia entre las energias del exciton y del 
modo fotonico en la microcavidad. La parte superior de la figura se hizo sin tomar en cuenta 
la interaccion de Coulomb (j3 = 0), de forma que los efectos de Coulomb son evidentes. Los 
dos picos en el panel inferior corresponden a resonancias entre el modo fotonico y la energia de 
uno u otro estado del par en el punto. 

2.4.3 Energia de excitones mult ip lenient e cargados a partir de la 
ecuacion de Bethe-Goldstone 

La ecuacion de Bethe-Goldstone es la base de la denominada aproximacion de pares indepen- 
dientes en nucleos [58J. Significa una mejora apreciable sobre los metodos de Hartree-Fock o 
de partfculas independientes en el sentido que se consideran las correlaciones de pares en la 
funcion de onda del micleo. Ha si do utilizada extensivamente en Ffsica Nuclear para el calculo 
de la energia de cohesion de la materia nuclear y de la energia (masa) de nucleos intermedios. 

En el trabajo [2U|, La ecuacion de Bethe-Goldstone fue reformulada para sistemas de excitones 
confmados en puntos cuanticos y empleada para el calculo de la energia de estos sistemas. Se 
investigaron sistemas desde el biexciton (2 electrones - 2 huecos) hasta conglomerados de 12 
pares, considerando ademas sistemas cargados, 4e - 2h y 6e - 2h. La diferencia E N+ i — E N 
entre las energias de N+l y N excitones permite evaluar como depende la frecuencia del laser 
que crea un nuevo par con el numero de pares presentes en el punto cuantico. De nuestros 
resultados se dedujo, ademas, la aparente inestabilidad del sistema 4e - 4h en un pozo cuantico. 

La ecuacion de Bethe-Goldstone describe el movimiento de un par de fermiones que ocupan los 
estados a y 7 por debajo del nivel de Fermi y se dispersan hacia estados por encima del nivel 
de Fermi. Dicha ecuacion se escribe: 

(Ti + T 2 + QaryVtyaj = E a ^ ai , (2.62) 
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Figura 2.12: Energias de sistemas de N excitones calculadas a partir de la ecuacion de Bethe- 
Goldstone. 

donde T son los operadores de energfa de partfculas libres en un campo externo y V es la 
interaction entre pares de particulas. Q ay realiza la proyeccion sobre los estados disponibles 

Q ai = \a 1 )(a 1 \+ A'Xa*', A'|. (2.63) 

Escribiendo: 

^7 = l«7)+ £ <VaV,A'>, (2-64) 

H',X>IX F 

obtenemos para C^J y E ai : 

= -(ti,\\V\a,>y), (2.65) 

£« 7 = eW+ e (°) + (a )7 |V|a, 7 ) 

+ £ (a, 7 |y|^,A')C;7 A ,. (2.66) 
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Las ecuaciones (|2.65j) constituyen un sistema lineal de donde se obtiene C®J como funcion de 
E ai , despues sustituyendo en (|2.fifij) llegamos a una ecuacion trascendente para Efxy. La energia 
total se halla de: 

En la Fig. l2.12l mostramos la energia de N excitones (2 < N < 12) calculadas por este metodo. 
El parametro f3 esta relacionado con la razon entre las energias caracteristicas de Coulomb y 
de oscilador, f3 = \J E c /(huj). Notese el escalamiento con respecto al mimero de pares y a la 
constante de interaccion. 

2.5 El metodo de Monte Carlo variacional en la descrip- 
cion de bosones en trampas y la proyeccion por Monte 
Carlo de la funcion BCS 

En los trabajos [Hj y [H] se utilizo el metodo de Monte Carlo para estimar, respectivamente, la 
energia de un sistema de alrededor de 200 bosones confinados en una trampa y para proyectar 
la funcion de onda BCS que se habia empleado en sistemas de excitones y en el calculo de la 
energia de apareamiento de micleos. 

2.5.1 Energia de sistemas de bosones cargados en trampas de atomos 

El perfil de densidad, el potencial qui'mico y otras magnitudes en sistemas de 10 4 - 10 6 atomos 
confinados en trampas y enfriados por laser hasta temperaturas muy por debajo de IK han 
sido medidas experimentalmente en fecha reciente jHH]. Inicialmente, el interes se centro en la 
realization experimental del fenomeno de condensation de Bose para atomos bosonicos (espin 
entero), aunque mas recientemente tambien se han estudiado los sistemas de fermiones alta- 
mente degenerados y los sistemas mixtos de Bose y Fermi. 

En el trabajo ^7j se calculo la energia del estado base de sistemas de hasta 210 bosones con 
interaccion de Coulomb entre ellos y confinados en una trampa cuasi-bidimensional. Aunque 
los sistemas de bosones con interaccion de Coulomb no han sido estudiados experimentalmente, 
la posibilidad de tener un plasma de particulas alfa en el micleo de una estrella neutronica [89J, 
o la posible existencia de bosones cargados en los superconductores a capas (HHI han motivado 
varios trabajos teoricos. El nuestro es el primero que estudio al sistema finite confinado en una 
trampa, aunque trabajos posteriores ya existen pT] . 

El calculo se hizo por dos vias: anahticamente utilizando los denominados aproximantes dobles 
de Pade de la Sec. I2.3l v a partir del metodo variacional de Monte Carlo (MC) jHZj. En realidad, 
la novedad del trabajo se inclina mas hacia la aplicacion de los aproximantes dobles de Pade 
a sistemas de cientos de bosones. En la presente seccion, sin embargo, nos concentramos en el 
metodo de Monte Carlo. 

El metodo variacional de MC combina el principio variacional de Ritz con la tecnica de MC 
para evaluar integrates multidimensionales. Siendo una funcion de prueba arbitraria, el 
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principio variacional nos dice que la energia del estado base satisface la desigualdad 



E gs < (9 T \H\* T ) 

< / d 2 n . . . & 2 r N \y T \ 2 E T (r X i ...,r N ) 



(2.61 



donde H es el hamiltoniano del problema y = ^j}H^t- La integral, de dimension 2N, 
es evaluada por MC. Para ello generamos puntos R = (fi, . . . ,r N ) con probabilidad dada por 
|\I/t| 2 - Tenemos: 



E, 



98 



' puntos 

El lado derecho de (|2.69J) es el estimado para E gs . 



(2.69) 
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Figura 2.13: Comparacion entre los resultados de Monte Carlo y los aproximantes dobles de 
Pade para 20 y 210 bosones. 
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El valor de N puntos utilizado en nuestros calculos fue 10 5 . La funcion de prueba contem'a 
correlaciones entre pares de particulas (factores de Jastrow) y sus tres parametros variacionales 
eran utilizados para minimizar la energfa. La generation de puntos en el espacio multidimen- 
sional se hizo a partir del algoritmo de Metropolis jJTj. Este es un algoritmo secuencial en el 
que el punto Rr+i se genera a partir de Rk de acuerdo con lo siguiente: 

(i) Se crea R' = Rk + SR, donde SR es una perturbation aleatoria. 

(ii) Se genera un numero aleatorio r distribuido uniformemente en (0,1). 

(iii) Si \^t(R')\ 2 /\^t(Rk)\ 2 > r, entonces Rk+i = R', en caso contrario Rr+i = Rk- 

Un ejemplo de los resultados de [Ej se muestra en la Fig. 12.131 El parametro (3 es funcion de la 
razon entre las energias caracteristicas de Coulomb y de oscilador: (3 3 = (me 4 / h 2 ) / (Huj) . El 
limite (3 —>■ corresponde a bosones no interactuantes, mientras que (3 — > oo corresponde a un 
regimen en que la repulsion de Coulomb domina y el sistema se comporta como un solido rigido 
(un cristal o molecula de Wigner). En j3 ~ 1 la contribution de las interacciones de Coulomb 
es altamente no perturbativa. La diferencia entre los estimados de Pade y por MC es inferior al 
1.5 % para 20 bosones y menor que el 4 % para 210 bosones en el rango del parametro mostrado 
en la Fig. l2~T51 



2.5.2 Proyeccion estocastica de la funcion BCS 

El metodo de MC en la variante de Metropolis fue utilizado tambien en el trabajo ^H] con 
el objetivo de proyectar la funcion de onda BCS al sector con un numero fijo de particulas, 
corrigiendo asi la mayor deficiencia de esta funcion cuando se aplica a sistemas flnitos, men- 
cionada en la Sec. 12.41 Esto no solo nos permitio hallar mejores valores para la energfa en el 
caso de sistemas de excitones sino que, de vuelta a la Fisica Nuclear, se calculo la energfa de 
apareamiento en diversos modelos de estructura nuclear y en el isotopo de Zr con masa atomica 
100. Los resultados nuestros son comparables en exactitud con los mejores calculos cuanticos 
por MC de esta magnitud 

La proyeccion de la funcion BCS sobre el sector de N pares se escribe 



n ^44 1 (2- 7 °) 

k=ji 



donde 



C 



N 



E 



J1, — ,JN 



-1/2 



Jn 



U 



31 



U 



(2.71) 



JN • 



El estimado para la energfa se obtiene de 



Ejbcs — (^bcs\H\^bcs) 

= J2 in. /; ./.vM./, ./v). (2.72) 

jl,-,3N 
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Figura 2.14: Energia de apareamiento entre niicleones en un modelo de dos niveles degenerados. 



donde e(j\, . . . , ]n) tiene interpretacion de energia cuando los estados ji, 
pados y los coeficientes 



W( Jl ,...,j N ) = C 



JV" 



3n 



U 



Ji 



u 



,]n est an ocu- 



;2.73) 



J AT 



que son explfcitamente mayores que cero, pueden interpretarse como factores de peso. La 
expresion ()2.72|) permite una evaluacion por MC, donde el conjunto de estados (ji, . . . , j^) se 
generan segiin Metropolis con peso W(ji, . . . , Jjv)- 

En la Fig. 12.141 se muestran algunos resultados en un modelo de estructura nuclear. En este 
modelo la energia de apareamiento se puede calcular exactamente, de modo que la razon entre 
la energia calculada y la exacta en el caso ideal es igual a uno. El grafico muestra esta razon 
como funcion de la constante de interaccion o apareamiento, G. La solucion BCS da un error 
de hasta un 40 %, mientras que LN da un error cercano al 15 % cuando G ~ 0.1. La proyeccion 
por MC a partir de la Ec. ()2.72|) . sin embargo, reduce el error relativo a menos del 5 % en todo 
el rango de variacion de G, llegando incluso a coincidir practicamente con la energia exacta 
cuando el apareamiento es fuerte. 
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2.6 La diagonalizacion exacta y el algoritmo de Lanczos. 
Aplicacion a puntos cuanticos pequenos 

En los trabajos [2111211123 I2E1 EI] se empleo la diagonalizacion exacta del hamiltoniano cuantico 
para obtener diversas propiedades de sistemas relativamente pequenos, con un maximo de 7 
particulas. 

En muchos de los problemas que se describen en esta y la proxima section, los calculos numericos 
tienen en comun el hecho de que aproximadamente el 80 % del tiempo computacional es dedi- 
cado al calculo de los elementos de matriz del potencial de Coulomb, (ni, n 2 |l/|r*i — r 2 | \na, n 4 ), 
donde los \n a ) representan estados del oscilador armonico en dos dimensiones. Para optimizar 
este tiempo, el autor y sus colaboradores calcularon todos los elementos de matriz no nulos 
entre estados de las primeras 20 capas de oscilador y los almacenaron en un archive Esto nos 
permitio, por ejemplo, reducir de varias horas a unos pocos minutos los calculos de HF para 
un sistema con 200 electrones. 

El primer paso en un proceso de diagonalizacion es la construction de la base de funciones 
multiparticulares. En nuestro caso, la base fue conformada por determinantes de Slater de N 
particulas: S = \ni, . . . ,n^). Como se menciono en la Introduction, la dimension del espacio 
de Hilbert crece exponencialmente con el numero de particulas. Imponiendo la conservation 
exacta de magnitudes como el momento angular y la proyeccion del espin total se logra reducir 
un poco la dimension de esta base. Y utilizando tecnicas como el algoritmo de Lanczos para 
hallar, al menos, un sector del espectro de energias pudimos abordar problemas con un maximo 
de 7 particulas. 

2.6.1 Absorcion y luminiscencia inter-banda en puntos cuanticos ba- 
jo campos magneticos muy intensos 

En el trabajo [21] se calculo el coeflciente de absorcion, la luminiscencia y el perfil de densidades 
de electrones y huecos en puntos cuanticos pequenos (de a 3 electrones) en presencia de 
campos magneticos muy intensos, entre 8 y 50 Teslas. Experimentos de luminiscencia en pozos 
cuanticos (donde existen muchos electrones confinados en una region cuasi-bidimensional) bajo 
campos magneticos tan altos han sido reportados recientemente [§3*] . 

Los estados finales de un proceso de absorcion e iniciales en uno de luminiscencia contienen 
un par electron - hueco adicional. El sistema mas grande considerado fue el 4e - lh. Las 
funciones de onda de los estados iniciales y finales fueron obtenidas por diagonalizacion exacta 
de la ecuacion de Schrodinger en bases de aproximadamente 5,000 funciones. Debido a que los 
campos son tan altos las energias convergen rapidamente aim cuando la dimension de la base 
es relativamente pequena. 

Las intensidades de las lineas de luminiscencia, por ejemplo, se calculan de la siguiente manera: 

/wHWW, (2-74) 

donde D = e nh>n es el operador dipolar interbanda. El estado n del hueco tiene momento 
angular invertido y proyeccion del espin invertida respecto del estado del electron. En el articulo 
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Figura 2.15: Intensidad del pico principal de luminiscencia correspondiente a excitones 
multiplemente cargados. 



se utiliza un modelo simple de dos bandas para puntos cuanticos de GaAs. 

Se obtuvieron resultados interesantes como la re-estructuracion de los estados de mas baja 
energia al variar el campo magnetico, dando como resultado el corrimiento de los picos de 
absorcion y de luminiscencia, asi como variaciones abruptas del coeficiente de absorcion o de 
la luminiscencia para valores de campos correspondientes a los llenados 1/3, 1/5, etc del efecto 
Hall cuantico. 

Como ejemplo de los resultados de este trabajo en la Fig. 12.151 mostramos los maximos de 
luminiscencia correspondientes a multiexcitones con un exceso de carga N e , es decir la recom- 
binacion desde sistemas con N e + 1 electrones y un hueco. Observamos un comportamiento a 
saltos de las intensidades cuando N e > 2. Los valores de campo en que esto ocurre corresponden 
a los llenados mencionados. 



2.6.2 Tunelamiento resonante a traves de un punto cuantico 

Experimentos recientes han mostrado la factibilidad de controlar el numero de electrones en un 
punto cuantico (desde uno hasta cuarenta), a la vez que reportan mediciones de la conductancia 
en los mismos [HI EH1 ESI E3 EH! • 

En el trabajo (21j se tomaron parametros que reprodujeran aproximadamente las condiciones de 
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los experimentos en j^H EI3 ESI EZI EH] y se calculo exactamente la conductancia y la densidad 
de niveles para un punto de 6 electrones. Hay que resaltar que el punto cuantico se halla 
separado de los electrodos por barreras de potencial, por lo que la conduction a traves del 
mismo se realiza por medio de un proceso de tunelamiento resonante. 

Las energias y funciones de onda de 6 y 7 electrones en el punto fueron obtenidas a partir de 
la diagonalizacion de la ecuacion de Schrodinger en bases de hasta 40,000 funciones utilizando 
el metodo de Lanczos para tratar matrices tan grandes. 
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Figura 2.16: Densidad de niveles de energia para un punto cuantico con 6 electrones en 
estados polarizados del espin. Las barras representan los resultados de la diagonalizacion exacta, 
mientras que los smibolos son estimados a partir del coeficiente de trasmision. 



El metodo de Lanczos hace uso del algoritmo de Gram-Schmidt para ortogonalizar una base de 
vectores que es creada por la propia matriz hamiltoniana (que se supone simetrica). Debido a 
que la base es creada por la propia matriz, las propiedades de convergencia ante truncamiento 
son muy buenas. Quiere decir, por ejemplo, que para obtener los primeros 30 autovalores con 
buena aproximacion basta representar la matriz en el subespacio de los primeros 500 vectores 
de Lanczos, aiin cuando la dimension de la matriz sea 40,000. 



Sea e*i un vector arbitrario. La primera iteration de Lanczos consiste en lo siguiente: 
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Calculamos ai = t\ - Hex, v 2 = Ht\ — aief, b 2 = \v 2 \ y definimos e 2 = v 2 /b 2 . Tenemos entonces: 



Hex = Oiei + b 2 e 2 . (2-75) 

En la segunda iteration calculamos: a 2 = e 2 ■ He 2 , V3 = He 2 — a 2 e 2 — b 2 ei, 63 = \vs\ y definimos 
63 = Tenemos: 



He 2 = b 2 ei + a 2 e 2 + b 3 e 3 . (2.76) 

Las siguientes iteraciones se hacen de manera analoga. Notese que en la base {e*j} el operador 
H se representa por una matriz tridiagonal. 

Nuestros calculos revelan que, a partir de una medicion precisa de la conductancia, se puede 
obtener la densidad de estados excitados en sistemas de pocos electrones. Esta afirmacion 
se ilustra en la Fig. 12.161 donde la densidad de estados obtenida por diagonalizacion exacta 
es comparada con la que se obtiene a partir de la conductancia (coeficiente de trasmision), 
calculada para un punto con 6 electrones y energias de excitation por debajo de 1 meV. En 
este intervalo, < AE < 1 meV, existen decenas de niveles de energia en el sistema estudiado. 
Nosotros consideramos solo los estados polarizados de espin, es decir, aquellos con espin total 
S z = N/2, que son los de menos energia cuando el campo magnetico es de alrededor de 10 
Teslas. 

En la figura se muestra un histograma con el mimero de niveles de energia en subintervalos de 
0.05 meV, el cual se compara con la magnitud: 



donde T peak es el valor del coeficiente de trasmision correspondiente a una energia de excitation 
AE y T^ es es el valor de este coeficiente cuando solo el primer estado excitado, que esta aislado 
del resto, esta en resonancia. Por ejemplo, cuando B = 12 Teslas el primer estado tiene AE m 
0.35 meV. La idea detras de la formula ()2.77|) es que al coeficiente de trasmision solo contribuyen 
los estados que estan en un intervalo de 0.05 meV alrededor de AE. El coeficiente de trasmision 
de cada canal independiente se puede estimar como T{ es y entonces el coeficiente total es 



T w 1 - (1 - T[ es ) n , (2.78 

de aqui la expresion (|2.77|) . 



2.6.3 Densidad de niveles de energia en un punto cuantico con 6 
electrones 

Motivados por los resultados de la section anterior, en el trabajo se estudio la densidad 
de niveles de energia en un punto cuantico con 6 electrones en un campo magnetico intenso. 
En estas condiciones, a energias de excitation por debajo de 1.5 meV y considerando todas las 
polarizaciones del espin, existen cientos de estados excitados. Notar que el rango de energias 
estudiados es pequeiio si se compara con la energia de ciclotron (ftw c ~ 12 meV), con la energia 
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Figura 2.17: Niveles mas bajos de energia del sistema de 6 electrones en cada sector con 
momento angular y espin total dados. 



caracteristica de Coulomb (E cou i = 8.4 meV) o con la energia de confinamiento del punto (ftw = 
3 meV). 

Se estudiaron todas las configuraciones de espin del sistema (desde S z = hasta 3) y se 
incluyeron estados de una particula de los primeros tres niveles de Landau para formar la base de 
funciones del sistema de 6 particulas. La matriz hamiltoniana en esta base alcanzo dimensiones 
de 850,000. Sus 200 autovalores mas bajos fueron obtenidos por medio del algoritmo de Lanczos, 
descrito mas arriba. 

En la Fig. 12 . 1 71 mostramos los niveles mas bajos de energia en cada torre con momento angular 
total dado, L. El campo magnetico se tomo igual a 7 T en el panel superior y 8 T en el inferior. 
Los saltos entre estados con S z = 3 y valores adyacentes de L son del orden de 0.6 meV, 
mientras que para cualquier otro valor de S z las variaciones son mas suaves. La contribucion 
de los niveles de Landau mas altos (segundo y tercero) a las energias es del orden de -0.4 meV, 
comparables con los saltos entre estados con L consecutivos. 

El resultado principal del articulo se muestra esquematicamente en la Fig. 12.181 El numero 
de estados en escala logaritmica es ploteado como funcion de la energia de excitacion. Se ve 
claramente que existen dos regimenes, por debajo y por encima de 0.3 meV. En cada sector, la 
grafica es muy bien aproximada por una dependencia: 



n(AE) = n Q exp 



AE 



(2.79) 



El parametro G cambia abruptamente en AE = 0.3 meV, senalando una discontinuidad en la 
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densidad de estados. Mediciones Raman bajo regimen de resonancia extrema pudieran poner 
de manifiesto esta discontinuidad. La discontinuidad en la densidad de estados se interpreta 
en terminos de la teoria de fermiones compuestos en el efecto Hall |102j . En los estados con 
AE < 0.3 meV, las cuasiparticulas forman clusters compactos donde la interaction entre ellas, 
a pesar de ser debil, debe considerarse. Por el contrario, los estados con AE > 0.3 meV pueden 
visualizarse como excitaciones del tipo Rydberg, donde uno o varios fermiones compuestos no 
interactuantes orbitan alrededor de un micleo de cuasiparticulas debilmente interactuantes. En 
la figura se ha incluido la prediction de esta teoria. 




i i i i i i i i i i i 

0,0 0,2 0,4 0,6 0,8 1,0 1,2 

Excitation energy [meV] 



Figura 2.18: Logaritmo del numero de niveles como funcion de la energia de excitation en un 
campo B = 7 T: los resultados exactos (cuadrados) se comparan con la ecuacion ()2.79|) (lineas 
discontinuas) y con la teoria de fermiones compuestos (lineas solidas). 



2.6.4 Transiciones opticas intra-bandas en el biexciton 

Las transiciones entre estados electronicos internos, causadas por absorcion de fotones en el 
infra-rojo lejano (5 - 10 meV), han sido estudiadas recientemente por medio de una tecnica, 
conocida como ODR, consistente en observar sus efectos en las lineas de luminiscencia [99|. 
Los experimentos fueron realizados en pozos cuanticos y la linea de luminiscencia tomada como 
referenda es la linea fundamental del exciton. 

En el trabajo se calculan los niveles de energia, funciones de onda y elementos matriciales 
del operador dipolar para el biexciton, defmido segun: 
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Figura 2.19: Energfas de excitation y probabilidades relativas de absorcion (solo las superiores 
al 10 % se representan) en el biexciton como funcion del campo magnetico. 



2 

D = J2(r ie -f th ). (2.80) 
i=i 

Campos magneticos entre y 5 Teslas son considerados. Los parametros utilizados (masas de 
las partfculas, etc) corresponden a los de los experimentos [99 . La mezcla de sub-bandas de 
Valencia se considero aproximadamente. El metodo de calculo fue la diagonalizacion directa 
de la ecuacion de Schrodinger en bases de hasta 10,000 funciones, para lo cual se hizo uso del 
algoritmo de Lanczos. 

Se hicieron predicciones para las posiciones e intensidades de las transiciones internas en el 
biexciton, las cuales deben ser observadas siguiendo la linea biexcitonica en la ODR. Se ob- 
tuvieron resultados interesantes sobre anti-cruzamiento de niveles y transferencia de fuerza de 
oscilador entre estados que colisionan. Un ejemplo de estos resultados se muestra en la Fig. 

em 



2.6.5 Absorcion intra-banda por excitones multiplemente cargados 



Consideremos un sistema de electrones cuyo movimiento esta confinado a un piano y sobre los 
cuales actua un campo magnetico perpendicular al piano de movimiento. Tambien a lo largo 
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de la normal, hacemos incidir luz polarizada circularmente. La absorcion de luz en el infrarrojo 
lejano mostrara un unico pico localizado en la denominada energfa de ciclotron, hio ce = heB/m e . 
Este resultado se conoce como Teorema de Kohn |l()()j y se debe a que la luz se acopla con el 
centro de masa del conjunto de electro nes. Para un sistema de huecos, que se comportan como 
partfculas con carga positiva, la posicion del pico de absorcion se halla por la misma formula 
(basta sustituir m e por m^), pero la polarizacion de la luz es opuesta. Los electrones absorben 
la polarizacion a + y los huecos la o~ . 

Este razonamiento permite comprender a grandes rasgos el experimento reportado en |101j . 
donde en esencia se mide la absorcion de ciclotron para electrones en un pozo cuantico como 
funcion de la concentracion de electrones en el mismo. A diferencia del problema mencionado 
en el parrafo anterior, un segundo laser con energfa superior a la brecha del semiconductor 
crea pares electron-hueco, con lo que el sistema no contiene solo cargas negativas, sino muchos 
electrones y un hueco. Por esa razon, la posicion del pico de absorcion deja de seguir la formula 
simple mencionada mas arriba y comienza a depender de la concentracion (densidad superficial) 
de electrones. Los resultados experimentales son "interpret ados" en terminos de excitaciones 
colectivas del tipo plasmon mas un hueco. Dichos resultados pueden resumirse en que el pico 
de absorcion se corre hacia el azul a medida que aumenta la concentracion de electrones. 

En el trabajo j3T] se calcula la absorcion en el infrarrojo de sistemas compuestos por unos pocos 
electrones (de 2 a 5) y un hueco. En el piano del movimiento se ha superpuesto, ademas, un 
potencial parabolico, con el cual modelaremos un punto cuantico |lUUj . El trabajo extiende los 
calculos presentados previamente para el biexciton (2 electrones mas 2 huecos) en un punto 

E2E|. 
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La comparacion directa de nuestros resultados con el experimento reportado en |l(Jlj desde 
luego que no es posible. El confinamiento lateral puede inducir comportamientos que no estan 
presentes en el caso infinito. Sin embargo, uno espera que cualitativamente se reproduzca el 
corrimiento al azul observado en el experimento. Por otro lado, no existen impedimentos para 
realizar mediciones directamente en un punto cuantico o en un arreglo de puntos. Notese que 
en [101J la absorcion en el infrarrojo se mide indirectamente a traves de los cambios inducidos 
en la luminiscencia cuando el haz infrarrojo es conectado. Y la luminiscencia es una tecnica 
extremadamente sensible y de facil implementation. 

Para los calculos utilizamos un modelo simple de dos bandas con parametros apropiados para 
el GaAs. La matrfz hamiltoniana se diagonaliza exactamente en una base de funciones de 
particular libres en un campo magnetico. Las dimensiones de esta matriz estan entre 40,000 
para los sistemas mas pequenos y 350,000 para los mas grandes. Matrices de estas dimensiones 
son diagonalizadas con ayuda del algoritmo de Lanczos, el cual nos permite obtener un conjunto 
de los autovalores mas bajos de energia. El error estimado para las energfas de excitation es 
de 0.02 meV y para las fortalezas de oscilador (normalizadas a la unidad) de 0.02. 

Un ejemplo de los resultados se muestra en la Fig. 12.201 En este caso, el sistema tratado es 
el trion o X~ ( es decir 2 electrones y un hueco, lo que daria una carga neta igual a -1. En 
general, X n ~ designara a un sistema compuesto por n + 1 electrones y un hueco). Los graficos 
estan separados de acuerdo a la proyeccion sobre el eje z del espm total de los electrones, S e y 
a la polarizacion de la luz absorbida. Los puntos en los graficos dan las posiciones de los picos 
principales de absorcion. Como information adicional damos tambien la fortaleza de oscilador 
en forma de un tirculo centrado en la posicion del pico de absorcion. De manera que un tirculo 
grande simboliza un pico fuerte de absorcion. 

En los graficos se ven tambien lmeas gruesas que simbolizan los maximos de Kohn para elec- 
trones o huecos en un punto cuantico parabolico 100 : 

AE^ = m e ± (2.81) 

AE { ± h) = hn hT ^, (2.82) 

donde Vt e = a/ uJq c + uS*J 4, u 0e es la frecuencia correspondiente al confinamiento lateral de los 
electrones, etc. El submdice ± en la energia se refiere a la absorcion de un foton con polarizacion 
a ± . 

El resultado principal de nuestros calculos se puede resumir en que la interaccion entre los 
electrones y el hueco provoca un corrimiento de los picos de absorcion con respecto a la posicion 
predicha por ()2.8H l2~%2|) . La dependencia de este corrimiento respecto a la carga neta, n, del 
punto cuantico es la siguiente: cuando n va de 1 a 3 el corrimiento disminuye (es decir, los picos 
se acercan a las posiciones dadas por (|2.81| l2~%2~|) ). mientras que para n = 4 crece nuevamente, 
manifest ando la tendencia observada en el experimento de que cuando la concentration de 
electrones aumenta los picos de absorcion se corren hacia el azul. 
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2.7 La aplicacion de los metodos de Hartree-Fock y de 
la fase aleatoria (RPA) a puntos cuanticos grandes. 
Resonancias dipolares gigantes y el efecto Raman 

Los trabajos [221 EI] estan dedicados a la absorcion intrabanda en sistemas de decenas o cientos 
de pares electron-hueco en un punto cuantico, mientras que los articulos [23 EDI E2J OH EH1 
38] abordan la dispersion inelastica de luz (efecto Raman) en puntos cuanticos. El comun 
denominador de todos ellos es que utilizan aproximaciones de campo medio tanto para el 
estado base (Hartree-Fock) como para los estados excitados (la denominada aproximacion de 
fase aleatoria, RPA). Mas abajo se describen algunos de los resultados mas interesantes de estos 
trabajos. 

2.7.1 Resonancias dipolares gigantes en sistemas neutros de elec- 
trones y huecos 

Las resonancias dipolares gigantes (GDR) son excitaciones colectivas de nucleos 103J que juegan 
un importante papel en la foto-absorcion y en otras reacciones nucleares. En el trabajo [22] 
se muestra que estados muy similares existen en sistemas de electrones y huecos en puntos 
cuanticos semiconductores. Cualitativamente, la GDR se puede ver como una oscilacion en 
contrafase de los electrones y los huecos. Es decir, los electrones moviendose en una direction 
y los huecos en la opuesta bajo la action del campo electrico oscilante de la luz. 

Las energfas de excitation y los elementos de matriz del operador de dipolo son calculados a 
partir de la denominada aproximacion de fase aleatoria (RPA) en sistemas de hasta 110 pares 
e - h, concluyendose que en las GDR se concentra aproximadamente el 98 % de la absorcion de 
luz con frecuencias de TeraHertz en estos sistemas. 

La RPA requiere como paso previo que el estado base sea obtenido en la aproximacion de 
Hartree - Fock (HF), por lo que se implemento la solution iterativa de las ecuaciones de HF: 

Et <E!»6« + P^ y E[(8,u\l/r\t,v) 

{ 7<[*F u,v 

- (s,u\l/r\v,t)]Cl u C^ (2.83) 

7<A4 U ' V ) 

donde los orbitales electronicos ocupados se expresan en terminos de las funciones de oscilador 
a traves de los coeficientes C ayS . Ecuaciones analogas a (j2.83j) se pueden escribir para la 
descomposicion de los orbitales correspondientes a los huecos. Note que los E^ son las energfas 
del oscilador en 2D 

Ef} = hyJu* + u*/4: {2k s + \l s \ + l} + ^l s 
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+g e fi B BSt, (2.84) 
E$ = ^-hy/^ + ^/A {2k s + \Q + 1} 

- — ^Is-gh^BSl (2.85) 
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Figura 2.21: Posicion de las resonancias gigantes y pigmeas como funcion del niimero de pares 
en el punto cuantico. 



Por otro lado, en la RPA el estado base \RPA) se supone que contiene correlaciones entre las 
particulas y los estados excitados se buscan en la forma 



V = Q^\RPA), 
donde el operador Q} viene dado por la expresion 



(2.86) 



Q ] = E^4e A + X h aX hlh x - Y^e\e a - Y x h a h[h a ). 



(2.87) 



El indice A corre sobre estados ocupados de HF, mientras que a corre sobre estados desocupados 
(por encima del nivel de Fermi). Los coeficientes X, Y se determinan por ecuaciones del tipo: 



EX A ee ~Y e j_ A e ^ _i_ R ee V e 

+ Bf x ^Xr}=^Kx, (2-88) 



CAPITULO 2. DESCRIPCION TECNICA ABREVIADA 



56 



en las cuales Kl son las energfas de excitation y las matrices A y B vienen dadas por 



Kir, = (E e<T -E eX )5 aT 5 Xfl + (3((a,fi\l/r\X,r) 

-((T,fi\l/r\T,\)) , 
A% Til = -i3{a,fi\l/r\X,r), 
B e A,,r = P((v,T\l/r\\,n)-{(T,T\l/r\ii,\)), 

B e X, T = -f3(a,r\l/r\\,n). (2.89) 

Usualmente, energias de excitation positivas (fisicas) y negativas (no fisicas) aparecen como 
soluciones de f!2.88|) . Las soluciones fisicas aniquilan el estado base 

Q\RPA) = 0, (2.90) 
y satisfacen la condition de normalization 

1 = E^I^aI 2 + \X h aX \ 2 - \Y x l\ 2 - \Y x h a \ 2 }. (2.91) 

Las ecuaciones ()2.88|) dan directamente las energias de excitation. Los elementos de matriz del 
operador dipolar, que causa la transition del estado base a los estados excitados, se expresan en 
terminos de los coeficientes X eY. La position de la GDR, es decir del estado que absorbe casi 
toda la fortaleza de la transition, como funcion del mimero de pares electron-hueco se muestra 
en la Fig. I2~2T1 

Por otro lado, en fecha reciente, se ha observado experimentalmente |104[ 1105^ 1106} IToT| H08j. 
que el 2 % de la foto-absorcion en nucleos se concentra en estados de mas baja energia de 
excitation. Por contraposition a las GDR, a estos estados se les ha llamado "resonancias 
pigmeas" . 

Los resultados del trabajo (221 muestran que la analogia con nucleos es casi perfecta pues los 
calculos de la absorcion en sistemas de electrones y huecos tambien dan como resultado una 
concentration del 2 % restante en estados de baja energia de excitation. (Ver Fig. 12.22}) . 

El numero de resonancias pigmeas, al igual que en nucleos, aumenta con el numero total de 
particulas. La position de la GDR y del promedio de las resonancias pigmeas tambien depende 
del numero iV de pares e - h en el punto cuantico. Esta dependencia es cualitativamente similar 
a la de nucleos, es decir la GDR se desplaza hacia energias mas altas al aumentar N, mientras 
que las resonancias pigmeas se desplazan hacia mas bajas energias. 

2.7.2 Dispersion inelastica de luz (Raman) en puntos cuanticos 

Consideremos, otra vez, el modelo simplificado de punto cuantico parabolico descrito en la 
section anterior, pero incluyamos en el decenas de electrones. La intention es acercarnos a las 
condiciones de los experimentos (Mj, en los cuales se miden las intensidades Raman en puntos 
con decenas o cientos de electrones. La dispersion de la luz en este modelo se puede calcular 
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Figura 2.22: Intensidades de las resonancias gigantes y pigmeas en puntos cuanticos. 



a partir de la teorfa de perturbaciones de segundo orden |lU9j . Es decir, el transito desde un 
estado inicial electronico, hasta un estado final, |/), con la consiguiente absorcion de un 
foton de frecuencia y emision de otro de frecuencia Uf, se realiza a traves del paso (virtual) 
por estados intermedios, \int). En las condiciones de los experimentos, la energia del foton 
incidente es hvi > E gap , donde E gap es la brecha del semiconductor, por lo que los estados 
intermedios contienen, ademas de los electrones iniciales, un par adicional electron - hueco. 



La amplitud del proceso se calcula a partir de la expresion: 

H^_ r \int, Ni — I) (int 
hvi — (E int — Ei) + iT int 



_ ^ (f,Nj- 1, l f \H+_ r \mt, Nj - I) {int, Nj - ljg^ji, N t ) 



hit 



donde H e _ r es el hamiltoniano de interaccion de los electrones con la radiacion y T int es el 
ancho energetico (fenomenologico) de los niveles intermedios. iVj es el numero de fotones en el 
haz incidente. 

El calculo de A^ requiere de : a) Las energias y funciones de onda uniparticulares de electrones 
y huecos, que son obtenidas en la aproximacion de Hartree-Fock. Las mismas son utilizadas 
como punto de partida en los pasos posteriores. La mezcla de sub-bandas de Valencia es 
considerada en un esquema de Kohn - Luttinger 4x4. b) Los estados finales de N e electrones, 
|/), hallados por medio de la aproximacion de fase aleatoria (RPA por sus siglas en ingles) [55] . 
c) Los estados intermedios de iV e + 1 electrones y un hueco, \int), que se obtienen a partir del 
denominado formalismo pp-RPA. Y, finalmente, d) calcular los elementos de matriz de H e _ r 
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con est as funciones y realizar la suma (|2.92|) . 

Con las amplitudes Aji uno calcula la section eficaz de dispersion: 



da 



~ 2_] \ A fi\ 2 $( E i + hlJ i -Ef-hu 



f 



fh 



(2.93) 



que es la magnitud medida experimentalmente. 

Como se menciono antes, en el conjunto de experimentos discutidos en jM] se obtienen los 
espectros Raman en puntos con decenas o cientos de electrones. Mediciones sin o en presencia 
de campos magneticos, sin o teniendo en cuenta la polarization de la luz incidente y reflejada, 
bajo diferentes condiciones de resonancia, son reportadas en estos experimentos. De forma 
simplificada, podemos resumir los resultados asf: bajo condiciones de resonancia extrema (es 
decir cuando hui practicamente coincide con E gap ) el espectro Raman es dominado por excita- 
ciones uniparticulares, mientras que a 40 meV o mas por encima de E gap los picos Raman estan 
asociados a estados finales que representan excitaciones colectivas (de espin o de carga). En los 
experimentos, hVi no va mucho mas alia de 40 meV por encima de E gap con el objetivo de no 
inducir procesos con fonones opticos. El efecto principal del campo magnetico es mezclar las 
excitaciones de espin y de carga, asi como desplazar la position de los picos Raman. 
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Figura 2.23: Intensidad del pico Raman asociado a un estado final monopolar de carga como 
funcion de Hva. Una resonancia en el canal de salida a 1642 meV es modelada. 



En el articulo [32] se presenta un analisis detallado de la dispersion Raman a cero campo 
magnetico. Se consideran tres condiciones de resonancia: a) E gap — 30 meV < hvi < E gap , b) 
E g ap < hvi < E gap + 30 meV (resonancia extrema) y c) E gap + 30 meV < hu^. El caso a) no ha 
sido abordado experimentalmente pero en el articulo se muestra su utilidad para identificar los 
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diferentes picos en el espectro. La teorfa presentada en [32] reproduce de forma cualitativa los 
aspectos mas relevantes observados en el experimento y senala otros, como por ejemplo el papel 
de T int y su dependencia con la energia de excitacion, que no habfan sido notados previamente. 

Como ejemplo y resumen de los resultados de en la Fig. l2.23l mostramos la intensidad del 
pico mas importante del espectro (excitacion monopolar colectiva de carga) como funcion de la 
energia de excitacion, hvi. La brecha efectiva en este caso es aproximadamente 1600 meV. En 
el caso hvi < E gap observamos una dependencia monotona de la intensidad del pico Raman. En 
resonancia extrema, por el contrario, las oscilaciones abruptas de la intensidad estan asociadas 
a estados intermedios resonantes con hvi. Cuando hui > E gap + 30 meV, el aumento de Y int 
hace que se pierda la dependencia oscilante con Y int observada antes y, en general, provoca una 
disminucion de la intensidad. Solo para determinados estados intermedios del tipo "exciton + 
excitacion colectiva", Y int conserva valores relativamente pequenos. En la figura, uno de estos 
estados es el responsable del aumento de intensidad en hvi « 1642 meV. En la figura se han 
incluido ademas, en forma de lmeas verticales, los numeradores que entran en la suma (|2.92)l . 
La comparacion con las intensidades Raman muestra que los efectos de interferencia son poco 
importantes en estos procesos. 

Otros resultados con y sin campo magnetico externo son presentados en [30] . El comportamien- 
to de las reglas de selection en presencia de un campo magnetico externo es estudiado en [38J. 
El caso de puntos cuanticos neutros o no dopados, donde una poblacion de electrones y huecos 
con N e = Nh es inducida por un segundo laser, es tratado en [2^] . Este ultimo caso hasta el 
momento no cuenta con una comprobacion experimental. 

2.8 La teoria degenerada de perturbaciones a segundo 
orden en puntos cuanticos de tamano intermedio 

La dimension de la matriz que representa al operador hamiltoniano crece exponencialmente 
con el numero de particulas. Asi, por ej., si consideramos un sistema de electrones en un punto 
cuantico cuasi-bidimensional y en presencia de un campo magnetico fuerte, donde la base de 
funciones mas idonea se construye a partir de los estados de un electron en un campo magnetico, 
sistemas con mas de 6 electrones se hacen practicamente inaccesibles a la diagonalizacion exacta. 

Consideremos, por ejemplo, 12 electrones con el espin polarizado a lo largo de B. Si tomamos 
78 estados en el primer nivel de Landau (ILL) para construir determinantes de Slater de 12 
electrones entonces, en el caso que la proyeccion del momento angular total es L = —132, el 
numero de funciones alcanza la cifra de 674 585. Esta cifra es manejable con el algoritmo de 
Lanczos segiin vimos en la Sec. 12.61 Sin embargo, intentemos incluir los efectos de los LLs 
mas altos. Estos efectos son importantes ya que pueden ser comparables en magnitud con 
las interacciones de intercambio. Tomamos 3 LLs con 78 estados en cada uno. El numero 
de funciones de la base es ahora mayor que 172 millones y la matriz hamiltoniana se hace 
inmanejable. 

En la presente section se muestra que el problema aiin se puede tratar utilizando la teoria 
degenerada de perturbaciones de segundo orden (PT2) |110j . Esto significa exactamente lo 
siguiente. El hamiltoniano se separa asi: H = Hq + V, donde Hq describe a N electrones en 
B y V — V con f + Vcoui contiene al potencial confinante y a la interaccion de Coulomb entre 
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Figura 2.24: Funciones de Slater en el subspacio de degeneracion, Sj, y en el orthogonal, Z, 
que entran en la suma de la Eq. ()2.95jl . El ejemplo es para un sistema con 4 electrones. 



electrones. 

Hq describe a un sistema cuantico altamente degenerado. Para el caso de 12 electrones en el 
ILL aparecen los 674 585 determinantes de Slater mencionados y a todos ellos corresponde la 
energia Eq = Nhu/2. A primer orden de teoria de perturbaciones degenerada (TP1) debemos 
diagonalizar la matriz 

H\f = (SilHlSA = EoSy + (SilVlSj), (2.94) 

donde los {Si} conforman el subespacio de degeneracion. En el siguiente orden perturbativo, 
TP2, debemos diagonalizar la matriz 

H*> = E a S,, + ( Si |l/ + Y, e! Z) eTz) ]S ' 1 <2J5) 

La suma sobre estados intermedios, {Z}, corre por el subespacio ortogonal, es decir, aquellos 
que cumplen (Z\Si) = 0. En los Z al menos uno de los electrones debe ocupar un estado de 
un LL mas alto que el primero. Vemos entonces que en TP2 la matriz a diagonalizar tiene la 
misma dimension que en TP1, pero el calculo de los elementos de matriz es mas complejo. 

Uno puede utilizar las propiedades de V para optimizar la suma sobre los Z. Por ejemplo, dado 
un Sj uno puede sumar solo por los Z tales que (Z\V\Sj) sea distinto de cero. Las propiedades 
de V a que nos referimos son: V con f es un operador de una particula, Vcoui es un operador 
de dos particulas y ambos conservan la proyeccion del momento angular total, L. Uno de los 
aspectos del trabajo consiste en disenar algoritmos que construyan estos estados intermedios. 
En general, la suma de la Ec. ()2.95|) se hara sobre estados de tres tipos: Z\, Z 2 y Z%. Un 
ejemplo para el caso N e = 4 que muestra las caracteristicas de estos estados se da en la Fig. 
IT2H 
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En el trabajo se reportan los resultados para el caso N e = 12. Los sistemas con N e = 2 (que 
se resuelve exactamente) y N e = 6 (cuyos autovalores de energia mas bajos fueron calculados 
bajo el esquema de Lanczos utilizando estados monoelectronicos de los tres primeros LLs 28 ) 
fueron utilizados como bancos de prueba para validar el algoritmo. La matriz tiene 10 
veces mas elementos no nulos que H^ l \ alcanzando tamanos del orden de 15 GB en disco duro. 
Su generation demora aproximadamente 2 semanas en un pequeiio cluster de 10 procesadores 
a 2.4 GHz. La diagonalizacion por el metodo de Lanczos toma de 1 a 2 semanas en un solo 
procesador. 

En la Fig. l2.25l mostramos los valores mas bajos de energia en cada sector con momento angular 
total L fijo. El campo magnetico se tomo como 10 T. Puede verse que la contribution de TP2 
a la energia es del orden de -0.6 meV, es decir aproximadamente -0.05 meV por electron. Estos 
valores son comparables a las energias de intercambio que marcan diferencias entre estados con 
distintos valores del espin total j2Hj. De aqui que si queremos determinar el espin total y el 
momento angular del estado base es necesario incluir los LLs mas altos en vez de reducir el 
calculo al ILL, como usualmente se hace |lllj . 
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Figura 2.25: Estimados por PT1 y PT2 de los niveles de energia mas bajos con espines 
polarizados de 12 electrones. B = 10 T. 



Capitulo 3 

Conclusiones y perspectivas 



Hemos presentado un conjunto de metodos no perturbativos y su utilization en problemas fisicos 
donde intervienen varias particulas. Estos problemas son dificiles de tratar y se distinguen 
del problema de una particula en un campo exterior, del cual existen numerosos ejemplos de 
soluciones analiticas, o de los problemas donde intervienen un numero muy grande de particulas 
y en los cuales son aplicables los metodos de la Fisica Estadistica. 

Como se menciono antes, los principales aportes nuestros no se concentran en los metodos, a los 
cuales hemos hecho modestas contribuciones, sino en el estudio de sistemas fisicos con varias 
particulas empleando estos metodos. A nuestro juicio, los resultados mas importantes de todos 
los resenados en la tesis son los siguientes: 

• El cuadro ofrecido por el metodo 1/D, que permite obtener expresiones analiticas para la 
energia del estado basico, energias de excitation y radios cuadraticos medios de sistemas 
que contienen desde 3 hasta alrededor de 10 particulas. En cierto sentido, la configuration 
rigida que se obtiene en la aproximacion principal del metodo es analoga a la orbita de 
Bohr en el atomo de Hidrogeno, mientras que las correciones provienen de relajar la 
rigidez y permitir oscilaciones. 

• La demostracion (numerica) de que algunos de los estados excitados en sistemas cuanticos 
de tres particulas (en especial, los denominados estados "respiratorios" o de "dilatation 
- contraction") son muy bien descritos como pequenas oscilaciones en un potential con- 
struido a partir de la funcion de onda del estado base (Sec. 12.2)1 . Algo similar es bien 
conocido por los especialistas que se dedican a la espectroscopia de nucleos atomicos y 
de otros sistemas: el espectro de excitaciones mas bajas muchas veces puede describirse 
como osciladores armonicos. 

• La construction de aproximaciones analiticas a la energia de sistemas con cientos de 
electrones basados en los aproximantes dobles de Pade. Se mostro que el error relativo 
de estos estimados no superaba el 2 - 3% cualesquiera fuera la fortaleza del potencial 
confinante. 

• La "demostracion", con varios hechos, de que los puntos cuanticos neutros fuertemente 
excitados, de forma que el numero de pares electron - hueco en ellos es relativamente 
grande, son, en gran medida, analogos a los nucleos atomicos. En efecto, en ellos existen 
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particulas de diversas especies, interacciones atractivas y repulsivas, la aproximacion de 
Hartree - Fock trabaja muy bien, la existencia de apareamiento entre fermiones, la for- 
mation de estados colectivos como las resonancias dipolares gigantes y otros hacen que 
la analogia sea casi perfecta. 

• Varios resultados especificos para puntos cuanticos con unos pocos electrones y huecos 
nos parecen muy interesantes y ameritan ser verificados en el experimento. Resaltemos, 
como ejemplo, el de la Sec. 12.6.21 referido al tunelamiento vertical a traves de un punto 
cuantico con 6 electrones acoplado debilmente a electrodos. En esa section se mostro 
que la medicion de conductancia permitiria estimar con buena exactitud la densidad de 
estados excitados en el punto. 

• Los calculos nuestros de dispersion Raman en puntos cuanticos con decenas de electrones, 
motivados por los experimentos realizados en los 90, los cuales no tienen similares en la 
literatura cientifica. 

• Y, fmalmente, la utilization combinada del metodo de Lanczos con la teoria de pertur- 
baciones a segundo orden, lo cual ha mostrado en un ejemplo especifico que se pueden 
tratar problemas cuya matriz hamiltoniana es de dimension extragrande. Esta combi- 
nation pudiera convertirse en una metodica de alcance mas general. 

Mirando hacia delante vemos que para nosotros existen varias perspectivas inmediatas de con- 
tinuar el trabajo. La primera consiste en introducir otros metodos en el analisis de sistemas 
con varias particulas. Entre los mas prometedores tenemos las integrates continuales evaluadas 
por Monte Carlo |112j . con las cuales el autor ya tuvo una primera experiencia en 1986 |113j . 
el metodo de los funcionales de densidad |114j . el denominado Grupo de Renormalizacion en 
Mecanica Cuantica j!15j y otros. 

La segunda posibilidad, para nosotros muy interesante, esta relacionada con el efecto Raman 
en puntos cuanticos multielectronicos. Ya mencionamos que en los ultimos anos los calculos 
de section Raman en puntos grandes son exclusivamente de nuestra autoria. Los experimentos 
habian practicamente desaparecido de la literatura desde el ano 98, pero en epoca reciente se 
nota un resurgimiento relacionado con la medicion del efecto Raman en puntos autoensam- 
blados |116t lllfj . Desde hace varios anos mantenemos una colaboracion con D.J. Lockwood 
en el Institute of Microstructural Sciences, en Ottawa, en el cual en estos momentos se estan 
realizando mediciones. La maquinaria nuestra de calculo esta lista, por lo que tenemos buenas 
posibilidades. Resaltemos que a esta tematica esta dedicada la tesis de doctorado en Fisica de 
Alain Delgado, proxima a defenderse. 

Tambien en la direction de vincular nuestros estimados teoricos con las mediciones experi- 
mentales esta la colaboracion que estamos comenzando con D. Whittaker en la University of 
Sheffield en el Reino Unido, que aborda la optica de puntos cuanticos acoplados resonante- 
mente a microcavidades semiconductoras. Herbert Vinck, cuya tesis de doctorado en Fisica 
esta en ejecucion, ha desarrollado una estancia larga en estos laboratories con resultados muy 
alentadores. Esperemos que la interaction con estos y otros grupos nos permita llevar nuestra 
comprension teorica de los problemas tratados a un nivel superior. 

Otra perspectiva de trabajo alentadora se refiere a la condensation de Bose - Einstein de 
excitones en puntos cuanticos. Los trabajos sobre condensation de excitones comenzaron en los 
anos 60, pero aiin no han recibido una confirmation experimental. Recientemente, la busqueda 



CAPITULO 3. CONCLUSIONES Y PERSPECTIVAS 



64 



se ha enfocado en la condensation en regiones finitas del espacio (puntos cuanticos) |118j . 
inspirandose en el exito que han tenido los trabajos experiment ales sobre condensation de 
atomos en trampas. En nuestro caso, ya esta implementado el algoritmo basado en la funcion 
BCS para tratar con decenas de pares confinados en un potential exterior. Faltaria ampliar un 
poco mas el tamaiio del sistema, de forma que los indicios de comportamiento macroscopico se 
hagan mas evidentes, y extender el formalismo al caso de temperatura no cero. El trabajo en 
esta direction ya ha comenzado. 

Finalmente, una extension natural de nuestro trabajo y que pudiera tener impacto consiste 
en abordar sistemas moleculares de interes biologico que involucren decenas o centenares de 
electrones. Los pasos que estamos dando en este problema van en dos direcciones: a) un 
tratamiento cualitativo que nos permita obtener estimados simples y rapidos y b) la imple- 
mentation de los metodos basados en los funcionales de densidad, los cuales representan en la 
actualidad la forma mas avanzada de abordar estos sistemas |120j . 
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